Mon. Not. R. Aston. Soc. 000,[T]— ?? (2012) Printed 10 December 2012 (MN MgC style file v2.2) 



ARGOS III: Stellar Populations in the Galactic Bulge of the Milky 

Way 

M. Ness 1 *, K. Freeman 1 , E. Athanassoula 2 , E. Wylie-de-Boer 1 , J. Bland-Hawthorn 3 , 
M. Asplund 1 , G.F. Lewis 3 , D. Yong 1 , R.R. Lane 4 , and L.L. Kiss 3 ' 5 ' 6 , 

^Research School of Astronomy & Astrophysics, Australian National University, Cotter Rd., Weston, ACT 2611, Australia 
2 Aix Marseille Universite, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, 13388, Marseille, France. 
'^Sydney Institute for Astronomy, University of Sydney, School of Physics A28, NSW 2006, Australia 
4 Departamento de Astronomia Universidad de Concepcion, Casilla 160 C, Concepcion, Chile 
5 Konkoly Observatory, MTA Research Centre for Astronomy and Earth Sciences, Budapest, Hungary 
6 ELTE Gothard-Lendiilet Research Group, H-9700, Szombathely, Hungary. 



Accepted 2012 TBD. Received 2012. As soon as possible; in original form 2012 February 



ABSTRACT 

We present the metallicity results from the ARGOS spectroscopic survey of the Galactic 
bulge. Our aim is to understand the formation of the Galactic bulge: did it form via mergers, 
as expected from ACDM theory, or from disk instabilities, as suggested by its boxy/peanut 
shape, or both? Our stars are mostly red clump giants, which have a well defined absolute mag- 
nitude from which distances can be determined. We have obtained spectra for 28,000 stars at 
a spectral resolution of R = 11,000. From these spectra, we have determined stellar param- 
eters and distances to an accuracy of < 1.5 kpc. The stars in the inner Galaxy span a large 
range in [Fe/H], -2.8 ^ [Fe/H] ^ +0.6. From the spatial distribution of the red clump stars 
as a function of [Fe/H] (Ness et al. 2012a), we propose that the stars with [Fe/H] > — 0.5 are 
part of the boxy/peanut bar/bulge. We associate the lower metallicity stars ([Fe/H] < —0.5) 
with the thick disk, which may be puffed up in the inner region, and with the inner regions 
of the metal-weak thick disk and inner halo. For the bulge stars with [Fe/H] > —0.5, we find 
two discrete populations; (i) stars with [Fe/H] w —0.25 which provide a roughly constant 
fraction of the stars in the latitude interval b = —5° to —10°, and (ii) a kinematically colder, 
more metal-rich population with mean [Fe/H] m +0.15 which is more prominent closer to 
the plane. The changing ratio of these components with latitude appears as a vertical abun- 
dance gradient of the bulge. We attribute both of these bulge components to instability-driven 
bar/bulge formation from the thin disk. We associate the thicker component with the stars of 
the early less metal-rich thin disk, and associate the more metal-rich population concentrated 
to the plane with the colder more metal-rich stars of the early thin disk, similar to the colder 
and younger more metal-rich stars seen in the thin disk in the solar neighborhood today. We do 
not exclude a weak underlying classical merger-generated bulge component, but see no ob- 
vious kinematic association of any of our bulge stars with such a classical bulge component. 
The clear spatial and kinematic separation of the two bulge populations (i) and (ii) makes it 
unlikely that any significant merger event could have affected the inner regions of the Galaxy 
since the time when the bulge-forming instabilities occurred. 
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1 INTRODUCTION 

Recent observations show that the stars in the Galactic bulge region 
comprise several kinematic and chemical components spanning 
a range in a— enhancement dBabusiaux et al.ll201Ct iBensbv et al .1 
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l2010h . Although the bulge region can no longer be described 
as a single stellar population, it appears to be a mostly old, 
a-enhanced system whose stars formed over a short period 
of time early in the life of the Galaxy dZoccali et all 120031: 
Melendez et al.ll2008l; rMcWilliam & Zoccalill2010l: llxcureur et al.l 



20071 : 1 Alves-Brito et alfcoiol) 



Near-infrared imaging l lOkuda et all dl977l) . COBE/DIRBE 
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jDwek et alj 1 19951 ; ISmith. Price. & Baker! |2004) revealed the 
boxy/p eanut morphology of the Milky Way bulge. Blitz & Spergell 
d 199 lh established the bar-like structure of the inner bulge region 
that is now understood to go with the boxy/peanut bulge struc- 
ture. The major axis of this bar points into the first Galactic quad- 
rant, at about 20 degrees to the Sun-cen ter line. Its semi-length 
is about 3.1 — 3.5 kpc i Gerhard 20021) and its axial ratios are 
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about 1: 0.33: 0.23 dDwek et alj|1995h. Near-IR sta r-counts (e.g 



lLopez-Corredoira. Cabrera-Lavers. & Gerhard 2005) provide evi- 
dence for a longer flatter c omponent of the bar, but the nature of 
this long bar is uncerta in. Marti nez- Valpuesta & Gerhar 
and iRomero-Gomez et al.1 1 20 111) argue that the apparent long bar 
is an artefact associated with leading spiral features at the end of 
the shorter primary bar. 

In C old Dark Matter (ACDM) simulations of Galaxy forma- 
tion (e.g. lAbadi et alj|2003l ; Ikobavashi & Nakasatdl201 it) , bulges 
are built up primarily through hierarchical mergers. However, it is 
now widely believed from simulations and observations that the 
boxy/peanut shaped bulges, like the bulge of the Milky Way, are not 
merger products but forme d via instabil i ty of t h e inner disk: see e.g. 
ICombes & Sanders! ( fl98lh ; lRaha et all dl99l!) ; lBureau & Freemar] 
1 1999!) . The early simulations showed that, after a few revolutions, 
the flat disk develops a bar which then becomes vertically un- 
stable and puffs up into a boxy/peanut-bulge structure. More re- 
cent N-body simula tions of bulge formation via disk instability 
d Athanas soulal 20081) help to understand the timescales and the spa- 
tial, kinematic and chemical abundance tracers of this instability 
process. For example, the simulations can be used to interpret the 
mapping of the early disk into the boxy/peanut-bulge via the insta- 
bilities. 

The simulations of bar/bulges produced by disk instability 
also show that cylindrical rotation is expected, and these are 
supported by observations of edge-on boxy/peanut bulges (e.g. 
iFalcon-Barroso et alj |2006). If the Galactic boxy/peanut bulge is 
a disk instability product, then we would expect it to show cylin- 
drical rotation. From their survey of M giants in the inner G alaxy 
and based on such kinematical criteria, [Howard et al. I J2009h con- 
strain any merger-generated component of the bulge to be less than 
8%. Several authors continue to argue in favour of the importance 
of dissipati onal collapse or mergers in the formation o f the Milky 
Way bulge feoccali et aUl2008l : lBlibusiaux et alj|2010l) . based on 
their interpretation of new chemical and kinematic data. 

We have completed a large spectroscopic survey of about 
28,000 stars in the bulge and surrounding disk in order to mea- 
sure kinematics and chemistry and so to constrain the formation 
process of the bulge. We aim to establish whether it is primarily 
a merger or primarily an instability product. We used models of 
classical bulges to estimate that 28,000 stars is sufficient to detect 
a 5% merger-generated bulge underlying an instability-generated 
bar/bulge. We note that this estimate is contingent upon the classi- 
cal bulge being a slowly rotating and dispersion-supported compo- 
nent. We would not be able to identify a classical bulge component 
which has been spun up into cylindrical rotation as in the simula- 
tions of : S aha, Mar tinez- Valpuesta. & Gerhard d2012h . 

This is the third paper in the ARGOS series and concentrates 
primarily on the metallicity distribution fu nction (MDF) for the 
ARGOS stars in the inner Galaxy. Paper I jNess et al]|2012al) ex- 
amined the nature of the bimodal distribution of apparent magni- 
tudes of the red clump stars along the minor axis of the bulge. This 
split clump appears to be a generic feature of boxy/peanut bulges 
and provides a useful guide to the interpretation of the metallicity 
components in the bulge. Paper II (Freeman et al., 2012) gave a de- 
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Figure 1. The 28 x 2° fields in our survey, in Galactic latitude and longi- 
tude. Filled circles indicate fields for which we have complete data. One 
field at (I, b) = (20°,— 7.5°) is incomplete; only 600 stars were ob- 
served for this field. The rectangle denotes the approximate extent of the 
boxy/peanut-bulge 



scription of the survey, our choice of fields, and the details of our 
observations, sample selection and analysis, including estimati on 
of stellar parameters and distances. Paper IV dNess et alj|2012cl in 
preparation) will describe the kinematical properties of the bulge 
as evaluated from the ARGOS stars. In Paper V, we will further 
interpret the origin of the MDF components by using our dynam- 
ical models to discuss how the instability of the disk generates a 
phase-space-dependent mapping of the disk into the bulge. 

In Section 2 of this paper, we examine the spatial and metal- 
licity distributions of stars in our sample, identify the components 
in the metallicity distribution in the inner Galaxy, and discuss the 
abundance gradients of these components in latitude and galacto- 
centric cylindrical radius Rq- In Section 3 we present the compo- 
nents that comprise our MDF and discuss their properties in Section 
4. We present our results for the a— enhancement of the bulge as a 
function of latitude in Section 5, and show the relationship between 
a— enhancement and [Fe/H] as a function of Rq - In Sections 6, 7 
and 8 we discuss the interpretation and origin of the MDF com- 
ponents that we have identified and compare our results to other 
studies. Section 9 gives our conclusions. 



2 RESULTS 
2.1 Overview 

For reference, Figure Q] shows the location of the 28 two-degree 
fields of the ARGOS survey. Each field has about 1000 stars, mostly 
red clump giants, as discussed in detail in Freeman et al. (2012). 

As part of the survey, we have obtained radial velocities for 
each star with an error < 0.9 km s _1 and determined the stellar 
parameters T cr r, log g, [Fe/H] and [a/Fe] for each of our stars. The 
[Fe/H] values are determined from the Fe lines in the Ca-triplet 
region of the spectrum, using a masked \ 2 technique against a 
grid of model spectra. We have excluded stars with poor quality 
or low signal-to-noise (S/N) spectra, giving a final sample of about 
25,500 stars. The resolution of the spectra is about 11,000 and the 
typical S/N for our stars is 75 per resolution element. From cali- 
bration against stars in well-observed globular and open clusters, 
the typical accuracies of the [Fe/H] and logg values are about 
0.09 and 0.30 dex respectively. Figure [2] shows an example of 
the (Teff, logg) distribution of stars for our minor axis fields a t 
(7,6) = (0°,-10°) with BaSTI isochrones dCassisi et al]|200fj|) 
of age 10 Gyr superimposed. The 6 panels show the distributions 
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Figure 2. (T cff , logg) values f or stars in the field (I, b) = (0°, -10°). 
overlaid with 10 Gyr isochrones {Cassisi et al. 2006). The number of stars is 
given in the left hand corner of each panel. The \ 2 fitting errors for log g 
are shown for each star and an error bar representing the calibration errors 
in log g along the y-axis and in temperature along the x-axis is shown in 
the right hand corner of the upper left panel. 



and isochrones for 6 intervals of [Fe/H] from —2.5 to +1.0. The 
isochrones have the mean abundance for the panel range in [Fe/H], 
and [a/Fe] = 0.4 for the isochrones at [Fe/H] < -0.25. The scat- 
ter of the stars about the isochrones is as expected. For some stars, 
the stellar parameters returned by our \ 2 process were too close to 
the edge of the model grid ( i.e. logg < 0.5) and these stars were 
excluded from the analysis. 

In order not to exclude metal-poor giants from the survey, 
we selected stars with (J — K) a > 0.38 and therefore expect 
to find some nearby dwarfs with spectral types later than about 
G2, as seen in Figure [2] These dwarfs are readily identified from 
their position relative to the isochrones in the (T c ff,log(?) plane. 
We see that most of the stars are clustered around the red clump 
(log. q= 2 to 3 and T e ff = 4500 to 5300 K) as we would expect 
i Zhao. Oiu. & Maoll200ll;IValentini & Munarill20ldlSaguner et all 



1201 lh . IZhao. Oiu & Mad ( 1200 ll) studied Hipparcos clump giants 
at high resolution and found two populations; a metal rich clump 
for which the gravity is higher ([Fe/H] > 0, log g= 2.5 to 3.1) plus 
a few more metal-poor stars with lower gravity ([Fe/H] < —0.3, 
log g = 1.9 to 2.4). We also find that the gravity of the main clump 
of stars is shifted towards lower values for lower [Fe/H], and we 
expect that the clump is not contributing significantly to the stellar 
population for [Fe/H] < —1.0. 

The clustering around the red clump region is seen across all 
of the fields in our survey. Figure[3]shows the density distribution 
in the (T c ff , log g) plane for the 25,500 stars in our survey, plotted 
in bins of metallicity against the isochrones. While the more metal- 
poor stars lie close to the isochrones, the more metal-rich stars de- 
viate progressively more and more from the isochrones. The devi- 
ations are larger for the giants that lie above the red clump, and are 
most marked for stars in the most metal-rich bin, [Fe/H] > 0.25. 
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Figure 3. Density images in (T fr , log g) space of 25 ,500 stars in our survey 
overlaid with 10 Gyr isochrones jCassisi et al.l2006l) . The number of stars is 
given in the left hand corner of each panel. The colorbars show the number 
of stars in each bin. 



We need to discuss these observed deviations from the isochrones, 
because we will be deriving isochrone distances for all of the sur- 
vey stars. The apparent deviations of our stars to the left may be due 
to several potential sources of error. Before discussing these error 
sources, we briefly summarise the relevant details of ou r procedure 
for deriving stellar parameters (see lFreeman et alj2012l for further 
discussion). 

While our \ 2 technique is in principle able to estimate T e e, 
log g and [Fe/H] from the AAOmega spectra in the Ca-triplet re- 
gion, we find in practice that there are too many degrees of freedom 
at our resolution, which leads to larger errors in the stellar param- 
eters. It is better to estimate r e g independently. This is co nsistent 
with the findings of lKirbv. Guhathakurta. & Sneder] d2008l) at their 
slig htly lower resolution of 8000. After c orrection for the redden- 
ing dSchlegel. Finkbeiner. & Davislll99^) , stellar temperatures are 
derive d from (J — K) a using a meta llicity-dependent calibration 
due to lBessell. Castelli. & Pled jl99^fl The procedure is iterative: 
starting with an initial estimate of [Fe/H] to get an initial estimate 
of temperature, new estimates of gravity and [Fe/H] are then deter- 
mined from the \ 2 technique against a grid of models. Typically 
four cycles are needed for convergence. The mismatch between 
isochrones and stellar parameters seen in Figure [3] for the more 
metal-rich stars could possibly be due to problems with the stel- 
lar models, the adopted reddening, or our temperature calibration. 
It seems unlikely that our abundance scale or our log g scale are the 
problem: they have both been carefully checked from observations 
of stars in previously well-observed star clusters covering a wide 
range of [Fe/H] and log g. 

The Ba STI isochrone s appear to work well for high metal- 
licity giants dCassisilfeoifjl) . Although we note that there is some 
evidence for a sprea d in age for the metal-rich dwarfs in the bulge 
dBensbv et aT]|2010l) . which could lead to some displacement to- 
wards higher temperatures for the corresponding metal-rich giants, 
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0.25 < [Fe/II] 1.1.) 




(a) (b) 

Figure 4. The most metal-rich stars on the isochrone from ICassisi et alj 
120061) . showing the effect of adopting a lower interstellar reddening (a) 
and an alpha-enhanced temperature calibration (b). 



we will assume that the 10 Gyr BaSTI isochrones are appropriate 
for our bulge stars. 

Some authors argue that the Schlegel red dening may be too 
high by 30 to 50% for regions of Ay > .5 jArce & Goodman] 
1 19991 : ICambresv. Jarrett. & BeichmarJl2005h . Although our fields 
were selected to have relatively low reddening, the regions of high- 
est reddening in our survey which are closest to the Galactic plane 
also correspond to the locations of the most metal-rich stars. On 
the other hand, (1) our stars all lie at Galactic latitude |6| > 4°; in 
the inner Galaxy. These stars are at least 300 pc above the Galactic 
plane, so the Schlegel reddening should be appropriate; (2) when 
selecting fields, we found that the Schlegel reddening was in excel- 
lent agreement with the reddening inferred from the morphology 
of the near-IR colour-magnitude diagram; (3) our most highly red- 
dened field at (I, b) = (0°, -5°) has a mean A v ~ 2.2. If the 
Schlegel reddening is systematically too high, then we would ex- 
pect to see a corresponding shift relative to the isochrones across all 
metallicity bins in this field. This was not seen. In summary, we be- 
lieve that systematic overestimation of the reddening is not respon- 
sible for the offset seen in Figure [3] at high metallicities. However, 
although a systematically high reddening is unlikely to be respon- 
sible for the offset seen in Figure[3]at high metallicities, an overes- 
timation of the reddening for the highest reddening fields only may 
explain the offset. Figure [4] (a) shows the position of the stars on 
the isochrone for the most metal-rich stars only ([Fe/H] > 0.25), 
where the temperature has been calculated for a reddening of 0.9 x 
E(J-K) only for E(J-K) > 0.3, where E(J-K) denotes the Schlegel 
reddening. The metal-rich stars now lie closer to the isochrone and 
the other metallicity bins are not affected. The effect on the derived 
MDF of using these lower temperatures is that the number of metal 
rich stars with [Fe/H] > would decrease by about 15% and the 
mean [Fe/H] of the MDF would decrease by 0.06 dex. 

We performed o ur analysis prior to the availability of the 
iGonzalez et al.l d201 lh redd ening map, which prov ides only partial 
coverage of our fields. The IGonzalez et alj d201ll) reddening val- 
ues in our fields are lower than the Schlegel values. For our field, 
(I, b) = (0°, —5°), the Gonzalez maps are about E(J-K) = 0.1 on 
average lower in estimated extinction compared to the Schlegel 
maps. This reduces the determined temperature by about 300 K 
(which corresponds to a lower mean [Fe/H] by about 0.2 dex). For 
the highest metallicity bin (Fe/H > 0), this lower photometric tem- 
perature gives a better fit to the stars on the isochrone, but at all 
other metallicities (Fe/H < 0), the stars are shifted too low in tem- 
perature and off the isochrone. 



The temperature calibration is the other possible source of the 
offset. Our adopted temperatures have been used for the calibra- 
tion of our parameters against a number of well-measured stars 
and clusters (see ARGOS II) and this temperature calibration is 
able to return the log g and [Fe/H] of our comparison stars within 
aiogg = 0.3 and oy Fe / H ^ — 0.09, respectively. The temperature 
relationships selected were for alpha-enhanced models for [Fe/H] 
^ -0.25 and for solar scaled models for [Fe/H] > -0.25. Over the 
T e ff range 4000 to 6000 K, a change of 50 K in T e g corresponds 
to a change of 0.05 in [Fe/H] and 0.15 in log g. A systematic tem- 
perature error as large as the offset seen in Figure [3] seems to be 
excluded (unless there is a systematic problem with the grid of 
model spectra such that the model spectra give the right [Fe/H] and 
log g values for stars whose temperatures have been overestimated 
from the phot ometry). We note however that, over our tempera- 
ture range, the lAlonso. Arribas. & Martinez- Roged dl999h temper- 
atures derived from (J — K) are typically up to 180 K cooler than 
the scale which we have adopted; the difference is largest for the 
more metal-poor stars. Figure [4] (b) shows the improved fit to the 
isoch rone by adopting the alpha-enha nced temperature relationship 
from iBessell. Castelli. & Pled d 19981) also for the metal-rich stars 
([Fe/H] > -0.25). The cooler temperatures move the metal-rich 
stars closer to the isochrone. The effect on the derived MDF us- 
ing these lower temperatures is that the total number of stars with 
[Fe/H] > decreases by 15% and the mean [Fe/H] of the MDF 
decreases by 0.03 dex. In summary, although the adopted temper- 
atures give accurate [Fe/H] and log g values, the temperature scale 
remains a possible reason for the offsets relative to the isochrones. 
Note that for the T c g - (J-K) relation, the J-K changes relatively 
slowly with T e g , so it is at more risk of errors than an alternative 
colour index such as V-K. However, our stars are sourced from the 
2MASS catalogue for which J-K is available and this colour rela- 
tion is also less sensitive to differential interstellar reddening which 
may be a problem for a few of our b = — 5° fields. 

Following this discussion, we formulate the following proce- 
dure for estimating stellar distances (see also paper II). For stars 
that are located near the clump (log g = 1.8 to 3.2 and T e ff= 4500 
to 53 00 K), we adopt an absolute magnitude Mk = —1.61 dAlvesI 
|2000|) which is appropriate for stars with [Fe/H] > —1.0 (95% of 
our sample). See the Appendix for further discussion of the abso- 
lute magnitude of the red clump adopted for distance calculations. 
As can be seen from Figure [2] the majority of stars are taken to be 
clump stars. For stars outside this range of T c s and log g, we derive 
Mk values using the BaSTI isochrone of age 10 Gyr and appropri- 
ate [Fe/H], adjusting only the temperature of the star to place it on 
the closest branch of the isochrone. 

For the red clump stars there are three main sources of dis- 
tance uncertainty: (i) the Schlegel reddening uncertainty, (ii) the 
2MASS uncertainty in the K magnitude and (iii) the red clump 
intrinsic magnitude dispersion. The re ddening uncertainty of the 
Schlegel maps is estimated to be 16% dSchelgel et al.lfl998T) . The 
mean reddening in our fields is E(B~V) = 0.2. A 15% reddening 
error would typically contribute a ok s = 0.02 to the K magnitude. 
In our stellar selection, we selected only those stars with 2MASS 
photometric errors 0.06 in each band. The mean error in each 
magnitude for the stars in our sample is much smaller, about 0.03 
mag. 

The red clump magnitude dispersion estimated by lAlvej 

J2000h is ok A = 0.22. Although this value includes the intrin- 
sic error of t he parallax measurement which has a om k = 0.1 
dAlvesll2OO0l) . we adopt the intrinsic dispersion of the clump to 
be ok ~ 0.22 to include the small contribution from reddening 
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and magnitude errors described above. As discussed in paper II, in 
the (logg Tee) plane, the red clump giants lie on a background 
of first ascent giants. The first giants with log g values close to the 
red clump have Mk values close to that of the clump. However, 
because our measuring errors for logg are relatively large (0.30 
dex), we have taken clump stars to be those with logg between 

I. 8 and 3.2, and our sample of clump giants therefore includes 
some brighter and fainter stars from the first giant branch (see pa- 
per II, Figure 7). For an unbiased population, most of these first 
giant branch stars are fainter than the clump stars. Our sample is 
not unbiased, because the 2MASS catalog is becoming incomplete 
for our fainter stars, and the numbers of first giant branch stars that 
we have included as clump stars is therefore reduced relative to the 
unbiased sample. 

For stars not on the clump, the distances are determined from 
the temperature and gravity via the the BaSTI l lCassisi et al]|2006l) 
isochrones. To determine its absolute magnitude, the star is shifted 
across in temperature to fit the nearest branch of isochrone. The 
error on log g (± 0.30) dominates the error of the absolute mag- 
nitude. For the 10 Gyr BaSTI isochrones, the absolute K magni- 
tude increases almost linearly with log g, with dA7d(log g) = 2.35. 
Therefore, <tm k = 2.35 x 0.30 = 0.70 for stars not on the clump. 
This corresponds to a distance error of about 38%. 

By using the red clump membership for as many of our po- 
tential clump stars as possible, across 1.8 < \ogg < 3.2 we are 
able to take advantage of the precision of the red clump absolute 
magnitude in calculating distances. Selecting stars over this inter- 
val has the advantage of including 99% of our red clump stars given 
our error of 0.3 dex in log g. The compromise in doing this is the 
contamination of non-clump stars that will be intrinsically brighter 
or fainter which will be bought into the inner region. For exam- 
ple, stars with a true log g at the upper limit of this selection will 
in reality lie closer to us than the bulge region which we define in 
our analysis to be Rq< 3.5 kpc. These stars will be moved fur- 
ther in distance from the Sun into the bulge. Similarly stars with 
a real \ogg at the low end of this selection will in reality be at 
larger distances than the bulge region. This kind of contamination 
by first ascent giants will be a prob lem for any distance e stimates 
using the red clump magnitude (e.g. Rattenburv e t al.l2007h . unless 
gravities are accurately know dValentini & MunarfeoiOh . Figure [5] 
shows the distribution of the stellar distance to the Sun as a function 
of log g for stars taken with distances from the isochrone only (at 
bottom) and using the red clump magnitude across 1.8 < log g < 
3.2 (at top) for our minor axis fields. These figures demonstrate how 
taking a range in log g as assigned clump stars to allow for our ob- 
servational errors, moves a fraction of closer and farther stars into 
the bulge region. The boxy structure seen in Figure[5]across xy re- 
flects the limits across our log g range that we have assumed for the 
clump and the K-magnitude range selected for our observations of 

I I. 5 to 14. The extent of the spread introduced in taking a range for 
the clump along the j/-axis is set by this magnitude range. We can 
reduce only our interval in log g, along the a>axis. We could have 
used a narrower interval in log g of say 2.0 - 3.0 but then we would 
have lost about 13% of the clump stars. Figure[6]evaluates the effect 
of our selection on the MDF for stars within Rq < 3.5 kpc on the 
minor axis using this interval in log g for red clump membership. 
The MDF at left is for the population where all distances have been 
taken directly from the isochrone. There are about 10,000 stars in 
the inner region using isochrone distances only. The centre panel 
shows the MDF generated for stars calculated with isochrone + red 
clump distances and there are about 14, 150 stars in the inner region 
using isochrone + red clump distances. The MDFs are only weakly 
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Figure 5. The distance, log g distributions for the ARGOS stars using dis- 
tances calculated for distances calculated from the isochrone + red clump in 
the interval 1.8 < log g < 3.2 (top) and from the isochrone only (bottom). 
Stars are binned into 100 equal width bins for each figure. 
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Figure 6. The MDFs for all stars on the minor axis within i?G< 3.5 kpc 
generated with the isochrone distances only (at left) and with isochrone + 
red clump distances (centre). The panel on the far right shows the fractional 
difference in the MDF due to the assumption of red clump magnitudes over 
this range in 1.8 < log g < 3.2 



affected. The far right panel shows the fractional difference of these 
distributions of MDF (isochrone + red clump) - MDF (isochrone 
only). The selection difference is small in the inner 3.5 kpc, and 
taking selections of smaller radii increases the relative number of 
metal rich stars and decreases the number of metal poor stars. For a 
selection cut of Rq < 1.5 kpc there are about 10% fewer stars in to- 
tal with -0.2 > [Fe/H] > -1.0 and 10% more stars with -0.2 > 
[Fe/H] > +0.5. 

To summarise, we take all of our potential clump stars with 
logg between 1.8 and 3.2 to have Mk = —1.61 ± 0.22 (about 
900 pc at the distance of the Galactic center), recognising that for 
some of the underlying first ascent giants the uncertainty in Mk 
will be somewhat larger. For the non-clump giants the distance er- 
ror is 0.70 mag, dominated by the uncertainty in the log g estimates. 
In our data analysis we verified that similar results were obtained 
for the smaller sample selection obtained using isochrone-only dis- 
tances. In Section 3 we evaluate Gaussian component fits to our 
MDFs. The relative fractions of the components are only slightly 
changed if we use isochrone-only distances. 

To evaluate the properties of the stars which lie in the inner 
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Table 1. Metallicity distribution of the 14,150 stars with |i? G | sg 3.5 kpc. 



Number of stars [Fe/H] range 



16 stars 


[Fe/H] sC -2.0 








84 stars 


-2.0 < [Fe/H] 






1.5 


522 stars 


-1.5 ^ [Fe/H] 


< 




1.0 


4219 stars 


-1.0 < [Fe/H] 






0.5 


6914 stars 


-0.5 < [Fe/H] 









2392 stars 


[Fe/H] > 









a? -1.0 
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all [Fe/H]: 100% 

[Fe/H] > 0: 17% 

> [Fe/H] > -0.5: 47% 
--- -0.5 > [Fe/H] > -1.0: 29% 

-1.0 > [Fe/H] > -1.5: 4% 
- - - [Fe/H] < -1.5: 0.91% 
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Figure 7. Radial distribution of ARGOS stars: In each panel, the x axis is 
the cylindrical distance Rq from the Galactic centre. Negative radii are for 
stars on the far side of the bulge and the vertical line represents the loca- 
tion of the Galactic center. Each star is plotted in the upper panel and the 
generalised and normalised histogram of the radial distribution is shown in 
the lower panel. The number of stars as a percentage in each [Fe/H] bin 
is indicated in the legend. The histograms are shown for all stars and for 
stars in bins of metallicity. The colours in the upper panels correspond to 
the colour coding in the lower panels. The [Fe/H] values are colour-coded 
from red (metal-rich) to blue (metal-poor). Note the changing line of sight 
distribution with metallicity and the dip near Rq = for some of the metal- 
licity bins. The dip near Rq = comes from the distribution of survey field 
centres; only the three fields on the minor axis contribute to the distribution 
near Rq = 0. The kernel size for the generalised histograms is a = 1.5 kpc, 
the largest estimated distance error. 

Galaxy, we take stars with cylindrical Galactocentric |.Rg| ^ 3.5 
kpc to be part of the inner Galaxy. This sample includes about 
14, 150 of the original 25, 500 stars. The metallicity distribution 
of these inner stars is summarised in Table[T] They are a mixture of 
all of the populations present in the inner Galaxy: the metal-poor 
halo, the thick disk, the thin disk around the bulge, and the bulge 
itself. 

2.2 The line of sight distribution of giants in the inner Galaxy 

Our sample includes giants all along the lines of sight with dis- 
tances up to ~ 55 kpc from the Sun on the far side of the bulge. The 
observed distribution of stars over distance is affected by our se- 
lection criteria, observational strategy and the stellar metallicities. 
Figure UJ shows the distribution of radii Rq for our entire sample 
of 25,500 stars. This figure does not reflect the underlying density 
distribution. Negative radii are for stars on the far side of the bulge, 
and the Sun is located at Rq = 8 kpc. The observed stellar density 
distribution is shown as the smooth curves (generalised histograms 
with kernel a — 1.5 kpc), for (i) all stars and (ii) stars in intervals 
of [Fe/H] to compare the relative density distribution of different 





Number ot stars 


0-2 


5100 


2-4 


12200 


4 — O 


DljU 


6-8 


2100 


> 10 


900 


> 15 


450 


>25 


100 


>35 


30 


s:3.5 


14,150 (55% of total) 


SC4.5 


19,650 (75% of total) 


\y\ S5 3.5 


17,350 (65% of total) 



Table 2. The number of stars in different intervals of radius Rq (about 
25, 500 stars). In the last line, y is a Cartesian coordinate along the Sun- 
center line, with origin at the Galactic center. 



metallicity groups to the total distribution of stars. In total, 60% 
more stars are observed per kpc along the line of sight on the near 
side of the bulge than on the far side. However, the line-of-sight dis- 
tribution is metallicity dependent. The majority of observed giants 
are located about 4 to 1 kpc from the Sun , with the peak about 3 
kpc in front of the bulge. iRobin et al.ld2012h estimates a peak in the 
thin disk number density at ~ 2.3 kpc from the centre of the Galaxy 
which may contribute to the larger number of metal rich stars seen 
at low Galactic latitudes on the near side of the bulge. The more 
metal-poor stars are less biased towards the near side of the Galac- 
tic center, with only 15% fewer stars observed in this range on the 
far side of the bulge. We preferentially sample the more metal-poor 
stars at larger distances because they are bluer, brighter and located 
at higher galactic latitudes. 

Figure [7j shows a minimum or inflection in the stellar number 
density around Rq = 0. This is due to the location of our fields; 
stars with \Rq\ < 1.0 kpc come only from our three fields at lon- 
gitude I = 0° and at latitudes b between 4 and 11°. This central 
feature is much less prominent for the more metal-rich stars, indi- 
cating that the more metal-rich stars are more concentrated to the 
inner regions of the bulge. The upper panel of Figure [7] shows that 
most stars with \Rq\ < 0.75 kpc have [Fe/H] > —1.0: only two 
stars from our sample in this region have [Fe/H] < —1.0. We will 
return to these properties of the inner region when we examine the 
metallicity distribution function (MDF) as a function of radius Rq 
(Figures 12 and 13). 

Although the lower limit of metallicity for our stars depends 
on radius (higher for \Rg\ < 0.75 kpc), the upper limit of the 
[Fe/H] distribution is similar at about [Fe/H] = 0.4 for all radii 
up to about \Rq\ = 5 kpc, with only a few stars above this level. 
Table [2] shows the distribution of the number of stars (giants) with 
Rq for stars seen out in the halo on the far side of the bulge. 

The most metal-poor stars in our sample, with [Fe/H] < —2.0, 
are located in the inner Galaxy but not in the innermost region. Of 
the 27 stars with [Fe/H] < -2.0, 18 lie between \R G \ = 0.9 and 
4.5 kpc. In the inner region, our most metal-poor star has [Fe/H] 
= —2.6 and lies at Rq= —3.8 kpc on the far side of the bulge. Our 
most metal rich star with [Fe/H] = 0.68 is at Rq = 2.0 kpc on the 
near side of the bulge. Although we find a few stars at metallicities 
[Fe/H] ~ 0.6, which are among the highest [Fe/H] values for bulge 
stars reported in the literature, we note that our \ 2 abundance scale 
is externally calibrated only for 0.1 > [Fe/H] > —2.1. Outside the 
range, the abundance scale comes directly from the grid of model 
spectra. 
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Figure 8. Histograms of metallicity for our three fields along the minor 
axis of the bulge at b = -5°, -7.5° and -f0° and with R G sC 3.5 
kpc. For (I, b) = (0°, -5°) 670 stars (73% of the total in the field) lie 
in this R G range with median [Fe/H] = -0.12. For (/, b) = (0°, -7.5°) 
there are 690 stars (73% of the total) and the median [Fe/H] = -0.34. For 
(I, b) = (0°, -f0°) there are 650 stars (67% of the total) and the median 
[Fe/H] = -0.46. The arrows indicate the median of the total distribution. 
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Figure 9. Metallicity Distribution Function for all stars within Rq ^ 3.5 
kpc. The vertical axis is the fraction of stars in each [Fe/H] bin. 



We have not included the solar neighbourhood dwarfs in Fig- 
ure|7] In total, our sample has about 2000 dwarfs located between 
about 1 and 3 kpc from the Sun. They have a median [Fe/H] 
= —0.1. These stars are not relevant to the present discussion of 
the inner Galaxy, and are not included in the discussion here. 



2.3 The [Fe/H] distribution in the inner Galaxy 

From the mean of the abundance distributions in our three latitude 
zones along the min or axis, we confirm the metallicity gradient 
seen in the bulge by IZoccali et al.l d2008l) . Our Metallicity Distri- 
bution Functions (MDFs) in the three minor-axis fields include a 
total of about 2000 stars with \ Rc,\ ^ 3.5 kpc (see Figure[8j, which 
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Figure 10. Fields with common latitude binned together for generating gen- 
eralised histograms of the metallicity distributions 



is about 70% of the whole sample of stars observed on the minor 
axis. Figure [9] shows the MDF for all of our survey stars with \ Rq\ 
^3.5 kpc. This distribution shows the relative proportion of stars 
in each [Fe/H] bin; it includes about 14,150 stars with contributions 
from fields across longitudes between ±20° and all four latitudes 
of b = -5°, -7.5°, -10° andfe = +10°. 

Our inferred MDFs shown in Figu r e [8] ar e qualitatively sim- 
ilar to those obtained by IZoccali et al.l d2008t) from their survey 
of about 800 K-giants observed at high spectral resolving power 
(R ~ 45,000) in their minor axis fields at b = —4°,— 6° and 
— 12°. However, our selection includes a larger fraction of metal 
poor stars due to our signi ficantly bluer color cut in our stellar se- 
lection compared to that of lZoccali et aljj2008l) . Quantitatively, our 
mean [Fe/H] values are about 0.15 dex lower on average, with our 
mean [Fe/H] decreasing from -0. 16 at b = —5° to [Fe/H] = -0.5 at 
b = — 10° . We find a somewhat lower vertical gradient of the mean 
abundance along the minor axis: —0.45 dex/kpc compared to their 
value of —0.6 dex/kpc. 



3 THE COMPONENTS IDENTIFIED IN THE [Fe/H] 
DISTRIBUTION 

We now examine the structure of the MDFs in our three latitude 
zones. To eliminate the effect of binning we examine the gener- 
alised histograms for our data in our fields, where each star is rep- 
resented by a Gaussian distribution centered on its [Fe/H] value 
and with a standard deviation of 0.13 (we adopted a kernel that is 
marginally wider than our estimated measuring error). To increase 
the number of stars, we have binned our data across all of the bulge 
fields, out to I = ±15°, for each latitude; b = -5°, b = -7.5° 
and b = —10°. The fields in the three latitudes which are grouped 
together are shown in Figure [lOl For a cut in radius of | Rg | ^ 3.5 
kpc, there are 4200 stars in the lowest latitude fields (b = —5°), 
2000 stars in the intermediate latitude fields (b — -7.5°) and 4000 
stars in the highest latitude fields (b = —10°). It was immediately 
obvious that the MDFs showed some structure which was qualita- 
tively similar over our three latitude intervals. 

Figures[TT](a), (b) and (c) show the generalised histograms at 
these three latitudes for stars within the bulge region |7?g 1 ^ 3.5 
kpc. The histograms show some structure and we can decompose 
them into a few Gaussian components. For all of the three latitudes, 
we can fit almost the same four Gaussian components, which we 
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denote A,B,C, plus an additional more metal-poor component 
E in the higher latitude fields (E is almost negligible at b = —5° 
and comprises only a very few stars at all latitudes). Figures[T"2l(a). 
(b) and (c) show conventional histograms of the metallicity dis- 
tribution in the three latitude zones. The individual MDF compo- 
nents can be identified. The smooth curves are gaussians with the 
same mean [Fe/H] as for the components in the generalised his- 
tograms of Figure QT] but with standard deviations a rea i reduced 
to account for the broadening effect of the generalised histogram 
kernel (ok = 0.13): the a rea i values for the gaussians shown in 
Figure [121 are given by a^ eal = <Jq H — a\, where ogh are the 
standard deviations derived in the multi-gaussian fit to the gener- 
alised histograms. 

Table[3]summarises the parameters for the MDF fits, corrected 
for the standard deviation of the generalised histogram kernel. We 
find slightly decreasing mean abundances of components A, B and 
C with increasing height below the plane |^|q 

Gaussian fits to the components have been adopted here for 
simplicity. It may be more correct to use skewed d istributions for 
the components (e.g. iGilmore. Wvse & Jones|[l995l) . This may re- 
move the need for the extra Gaussian components at the low metal- 
licity tail of the distribution, for which there are only a few stars. 
Components D and particularly E in the metal-poor tail are fit 
here only to demonstrate the changing fraction of metal-poor stars. 
Where stars with Fe/H <-1.5 are present, component E is included 
to highlight their relative fraction, although this component is gen- 
erally not statistically significant. It is instructive however to see 
how this metal-poor contribution changes with latitude and longi- 
tude. Table[3]demonstrates how the fraction of metal-poor stars in- 
creases with latitude from an almost negligible value in our lowest 
latitude zone. 

A Monte Carlo analysis was performed to determine the error 
in the MDF component parameters due to the fitting and sampling. 
We took the observed sample of [Fe/H] values and generated a new 
value of each [Fe/H], drawn from a Gaussian distribution around 
the mean of the data point and its standard deviation (we adopted 
a a = 0.13, fractionally larger than the measuring error). Revised 
Gaussian fits were then made to the generalised histogram gener- 
ated from this new distribution and the process was repeated 10 3 
times to determine the variance in the mean, amplitude, and stan- 
dard deviation of the 5 components. Table [3] summarises the out- 
come: for the mean [Fe/H] of the three main components A,B,C 
shown in Figure QT](a)-(c), the errors are ss 0.01. For the ampli- 
tudes they are typically 0.01, and for the standard deviations of 
the Gaussian components, they are about 0.01 to 0.05. The corre- 
sponding Gaussian components along the minor axis fields only are 
shown in Table|4] 



3.1 Justification of components 

The maximum likelihood method was employed using t he Ex pec- 
tation Maximisation (EM) algorithm of iMeng & Rubinl i ll 9931) for 
multi Gaussian fits within the statistical package B0, to determine 
the log likelihood for N=2,3,4,5,6 components for our MDF 



Gaussian Components Fit to Data 



-5° 





Mean [Fe/H] 


Weighting 


^real 


A 


+0.12 ±0.01 


0.38 ±0.01 


0.12 ±0.01 


B 


-0.26 ±0.01 


0.40 ±0.01 


0.12 ±0.01 


C 


-0.66 ±0.02 


0.20 ±0.01 


0.15 ±0.05 


D 


-1.16 ±0.02 


0.01 ±0.01 


0.15 ±0.02 


E 


-1.73 ±0.08 


0.002 ±0.01 


0.10 ±0.04 


-7.5° 




Mean [Fe/H] 


Weighting 


^real 


A 


+0.11 ±0.01 


0.16 ±0.02 


0.10 ±0.02 


B 


-0.28 ±0.01 


0.50 ±0.02 


0.16 ±0.02 


C 


-0.70 ±0.02 


0.30 ±0.02 


0.14 ±0.02 


D 


-1.19 ±0.09 


0.04 ±0.01 


0.12 ±0.05 


E 


-1.68 ±0.3 


0.006 ±0.01 


0.06 ±0.10 


-10° 




Mean [Fe/H] 


Weighting 


^real 


A 


+0.08 ±0.02 


0.12 ±0.02 


0.11 ±0.02 


B 


-0.30 ±0.01 


0.44 ±0.02 


0.14 ±0.01 


C 


-0.69 ±0.01 


0.38 ±0.02 


0.15 ±0.02 


D 


-1.21 ±0.03 


0.05 ±0.01 


0.13 ±0.02 


E 


-1.70 ±0.10 


0.01 ±0.01 


0.15 ±0.04 



Table 3. The mean, a and weights of the rive metallicity components iden- 
tified in the Galactic bulge in the three latitudes out to I = ±15°, for \Rq 
< 3.5 kpc 



Gaussian Components Fit to Data 



-5° 





Mean [Fe/H] 


Weighting 


°Y e a I 


A 


+0.15 ±0.02 


0.42 ±0.02 


0.11 ±0.02 


B 


-0.25 ±0.02 


0.42 ±0.02 


0.15 ±0.02 


C 


-0.71 ±0.02 


0.14 ±0.02 


0.13 ±0.02 


D 


-1.18 ±0.1 


0.01 ±0.005 


0.08 ±0.05 


E 


-1.80 ±0.1 


0.003 ±0.002 


0.11 ±0.07 


-7.5° 




Mean [Fe/H] 


Weighting 


°Y e a I 


A 


+0.13 ±0.02 


0.19 ±0.02 


0.08 ±0.02 


B 


-0.26 ±0.03 


0.46 ±0.02 


0.18 ±0.02 


C 


-0.70 ±0.02 


0.31 ±0.02 


0.16 ±0.02 


D 


-1.19 ±0.1 


0.04 ±0.005 


0.16 ±0.05 


E 


-1.68 ±0.1 


0.003 ±0.002 


0.08 ±0.07 


-10° 




Mean [Fe/H] 


Weighting 




A 


+0.09 ±0.03 


0.12 ±0.03 


0.12 ±0.01 


B 


-0.31 ±0.02 


0.41 ±0.02 


0.15 ±0.02 


C 


-0.70 ±0.02 


0.40 ±0.02 


0.15 ±0.02 


D 


-1.17 ±0.03 


0.06 ±0.01 


0.11 ±0.02 


E 


-1.68 ±0.2 


0.01 ±0.01 


0.08 ±0.06 



2 see Section 3. 1 for justification of number of components 

3 By fitting all Gaussian parameters in the more general method used 
here, we find sl i ghtly decreased gradients relative to those reported in 
iNess & Freeman! [20121) 

4 http://www.r-project.org/ 



Table 4. The mean, a and weights of the five metallicity components iden- 
tified along the minor axis inside Rq ^3.5 kpc 
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[Fe/H] [Fe/H] [Fe/H] 

(a) I ±15°, 6 =-5° (b)«±15°,fe= -7.5° (c) Z ± 15°, b = -10° 

Figure 11. Panels (a),(b),(c) show the generalised histograms for stars in the inner Galaxy with Rq ^ 3.5 kpc at latitudes b = —5° , b = —7.5°, b = —10°. 
Our smoothing kernel is a = 0.13. The number of stars in each panel is 4200, 2000 and 4000 stars respectively. Our Gaussian decomposition of the MDFs is 
shown; A - E. 




distributions given our probability distribution function. 

The probability density f(x) that is most likely to have generated 
the data points is represented by a mixture of k = l....n Gaussian 
functions G k \ with parameters 6 k = (/j, k , u k ) and mixture weights 
Pk so the equation takes the form 



f(x;e) = J2(p k G k ) 



where each Gaussian function is: 

G k = G(x;fj, k ,a k ) = ^^e »\ °* ) 



2na k 



(1) 



(2) 



where the probability distribution integrates to 1 and the weights of 
the components sum to 1 . 



J X I 1 



(3) 



x is the data of length N (number of [Fe/H] data points) 

The maximum log likelihood is then defined as 

L(x;e) = Yl? =1 f(x i; e) (4) 

and in our case where we have a Gaussian mixture model 



N n 

L(x;0) = ]~[ y^^pkGjxj- fi k , a k ) 
The log likelihood becomes 

N n 

lnL(x;0) = y^Jny^p k G(xj- fi k ,a k ) 



(5) 



(6) 



The EM algorithm maximises the likelihood function (minimises 
the negative log likelihood) and so returns a single value given 
a Gaussian mixture distribution. To assess the relative fits of the 
models, the Bayesian Information Criterion (BIC) was then applied 
which compares the relative fit of the models by incorporating a 
penalty on the maximum likelihood: 



BIC = -2 \nL + n ln(JV) 



(7) 



where iV is the number of data points and n is the number of free 
parameters which are being estimated. 

The Bayesian Information Criterion (BIC) method of assessing 
goodness of fit, which uses the maximum likelihood (and associ- 
ated Akaike Information Criterion (AIC)) are simplified metrics 
from a fully Bayesian analysis (Heavens 201 1). These methods 
find only a maximum and not average likelihood, however. The 
complete Bayesian approach calculates the likelihood over the full 
parameter space, thus assessing the predictive power of the model, 
given a prior. However, for our purposes, the information theory ap- 
proach is sufficient to find the best model. This is because, although 
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we aim to demonstrate that our fits are valid, our decomposition is 
also motivated by the peaks we can see in the smoothed distribu- 
tion, the behaviour we find for the red clump stars along the minor 
axis dNess et al]2012al) and the clear dependency of kinematics on 
rFe/Hl dNess et alj|2012cl in preparation) 

The BIC test returns a minimum (optimal) value for a 4- 
component fit (A,B,C,D) at b = —5° and b — —7.5° and a mini- 
mum for a 5-component fit (A,B,C,D,E) at b = —10°. In the Fig- 
ures QT|[T3] [14] and Tables [3] and [4] the component E at the low- 
est [Fe/H] group is included to demonstrate the changing number 
of the most metal-poor stars seen in the sample as a function of 
latitude, even though this is not a statistically significant fit. The 
maximum likelihood fitting routine returns the parameters of the 
X 2 minimisation fitting of Table [3] (except at the metal-poor tail) 
for the MDFs at b — — 7.5 and b = —10°. For these mixture distri- 
butions, component D is returned with weighting w 0.01 and with 
a wider standard deviation than for the \ 2 fit (c ~ 0.5). This dis- 
crepancy does not affect the other Gaussians fits nor does it alter 
the conclusions from the decomposition. The agreement of the two 
approaches indicates that the majority of stars are in components 
A,B,and C and these components are securely fitted. The metal- 
poor components D and E involve relatively few stars and their 
contributions for the different fields can be seen in Table [3] They 
are not the main target of our discussion here. 



3.2 Consequence of mixture model 

The component fitting described above identifies the three ma- 
jor Gaussian components A,B,C which are directly visible in the 
MDFs, plus the two minor components D and E. We ar e motivated 
here no t only by these fits but also by the split clump (see lNess et al.l 
2012a) and kine matics which are c learly different for our identified 
components (see lNess et al.ll2012cl in preparation). Independent of 
our MDF component fits, when we examine the kinematics and 
distribution of the stars partitioned into [Fe/H] intervals of 0.5 dex 
from +0.5 to —1.5 (these intervals correspond approximately to 
our components A-D), these subgroups show different kinematics 
and different split clump characteristics. 

If we were to adopt two components for our data rather than 
our best fit representation of 4 or 5 components, the optimal two 
components comprise a narrow metal rich distribution and a wider 
metal poor distribution, with (jui,cri) = (0.15,0.1), and (/X2,C2) 
= (-0.3, 0.35). This is similar to the compon ents identified by 
Ifflll et al.l d201 lh and bv lBabusiaux et alj | |2010[) in their b = -4° 
field. Both argue for a two-component bulge population, compris- 
ing a metal-rich population and an old spheroid. The number of 
components one considers clearly plays a important role in the in- 
terpretation. 



components B and C still dominate and component A is now quite 
weak. The metal-poor components increase in significance with in- 
creasing | b |, and more metal-poor stars ([Fe/H] < — 1) appear at 
higher latitudes. 

In preparation for the discussion in Section 7 about the nature 
of the metallicity components seen in the inner Galaxy, we show 
MDFs for intervals of Rg which demonstrate the changing con- 
tribution fraction of the different components with distance from 
the Galactic center. Figure [T3l shows the MDF components in var- 
ious distance intervals as indicated below the individual panels for 
b = —5°, and Figure [T4lis similar for b = —10°. The data for each 
latitude zone are combined as before for I — ±15°. 

Figures[T3landll4lshow that the components A - E are present 
at both latitudes out to |.Rg| ~ 5 kpc. The widths of the individual 
components are approximately unchanged, but their contribution 
fraction changes with latitude and Rg. Examined separately, the 
near and far sides along the line of sight show comparable MDFs 
out \Rg\ < 4.5 kpc. There are fewer metal rich stars observed on 
the near side however and a radial change in the morphology of the 
MDF occurs for |-Rg| > 4.5 kpc: at these larger radii, components 
A and B merge (see for example components A and B in panels 
e and / of Figure [TSl . Note that this conclusion depends on the 
adopted metallicity gradient for component A as discussed below. 
The distance cuts show that within Rg < 4.5 kpc the two compo- 
nents A and B are necessary to fit to the stars with « [Fe/H] > 
-0.5. However, outside the bulge, (for a radius of Rg > 8kpc at 
b = —5°, and Rg > 4.5 kpc at b = —10°), a single component is 
sufficient to fit to the stars with [Fe/H] > -0.5 (see Figures [T3"l(i) 
andEKi)). 

At b = —5°, component B dominates between 0.75 to 3.0 
kpc but inside of this region the most metal-rich component A is 
marginally strongest. Out to 4.5 kpc, component C is no more than 
about 50% of that of A. At b = —10°, component B is the strongest 
out to \ Rg\ = 6 kpc and component C is a larger fraction of com- 
ponent B out to 4.5 kpc compared with the lower latitude. Within 
Rg < 4.5 kpc, component A is relatively weak, at no more than 
about 30% of the strength of C. 

At b — —10°, component A which appears as a small fraction 
in the inner region is no longer present on the far side of the bulge 
(see Figure [141(g) and (h). At both latitudes, component B appears 
in the largest fraction on the far side of the bulge and component A 
is the largest fraction on the near side of the bulge outside of |J2g| 
> 4.5 kpc. 

Components D and E include only a small number of 
stars. At & = —5°, their fraction increases with Galactic radius. 
The fraction of very metal-poor stars (in D and E) is higher at 
b = — 10° , particularly near the Galactic centre (| Rg | < 0.75 kpc). 



4 PROPERTIES OF THE COMPONENTS 

The decomposition of the MDF into the A,B,C (D and E) compo- 
nents allows us to quantify the varying contributions of the metal- 
licity components over I and b and also over distance along the line 
of sight. 

Figure QT] shows how the relative contributions of the com- 
ponents change with latitude. At & = —5°, the more metal-rich 
components A and B dominate, with a smaller contribution from 
component C and only a few stars from the metal-poor components 
D and E. At b = —7.5°, components B and C are the strongest; 
component A is clearly weaker than at b = — 5°. At b = —10°, 



4.1 Comparison with disk fields 

One of our fields lies out in the disk at (I, b) = (—31°, —5°), for 
which the minimum value of \Rg\ is about 4 kpc. The same com- 
ponents A,B,C are seen in this disk field (Figure \T5[. From Figures 
[T3]and[T4] these three components are also seen in the radial inter- 
vals extending out into the disk at smaller 1 1 |. The presence of the 
same three components A-C in the bulge and the inner disk makes 
sense. In the inner disk, these components are present in their undis- 
turbed form, as disk stars, while in the bulge fields the stars of these 
components have been incorporated dynamically into the bulge. 

With increasing distance from the Galactic center, we note 
again the change in component A relative to components B and C. 
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(a) \R G \ <0.75kpc 



(b) \R G \< 1.5 kpc 



(a) | Rq I < 0.75 kpc 



(b) |H G | < 1.5 kpc 





(c) 1.5 kpc < \ R G \ < 3 kpc 



(d) 3 kpc < | R G | < 4.5 kpc 



(c) 1.5 kpc < |ii G | < 3 kpc 



(d) 3 kpc < | R G | < 4.5 kpc 





(e) 4.5 kpc < | Rq I < 6kpc 



(f)5kpc< |8 G <7kpc 



(e) 4.5 kpc < | Rq I < 6kpc 



(f)5kpc< |R Gn „ r j <7kpc 





(g) 4.5 kpc < |R G , <6kpc 



(h) 5 kpc < | H G | < 7 kpc 



(g) 4.5 kpc < |H G , <6kpc 



(h) 5 kpc < | H G | < 7 kpc 





(i) 8 kpc < | H Gfar I < 10 kpc 

Figure 13. MDF decompositions for different radial intervals alb = —5°. 
The number of stars in panels a-i are 300, 1 100, 2300, 2000, 600, 200, 120 
and 100, 70, respectively. 



(i) 8kpc< |H Gfar I < 10 kpc 

Figure 14. MDF decompositions for different radial intervals at b = —10° . 
The number of stars in panels a-i are 180, 850, 2450, 1700, 250, 150, 200 
and 120 and 50, respectively. 
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(a) I: 



-31 c 



-5° 



Figure 15. MDF for the field (I, b) = (-31°, -5°) from 440 stars with 
\Rg\ - 4.0-4.9 kpc. The Gaussian decomposition of the MDFs is shown. 
We see the same components out in the disk as in the bulge fields 



For | Rq I > 4.5 kpc, components A and B merge and, at b — — 10° , 
four components (B, C, D and E) are sufficient to fit the data, as- 
suming a radial gradient in the mean abundance for A. At 6 = — 5° , 
where components A and B merge outside of the inner bulge re- 
gion, component A is dominant close to the Sun. On the far side at 
Rg > 4.5 kpc, component B is dominant and components A and C 
is roughly equal out to Rg = 7.0 kpc. At b = —10°, component A 
is similarly the most significant near the Sun and at this higher lat- 
itude on the far side of the bulge, for Rg > 4.5 kpc, component C 
is the strongest. Table[5]shows the relative heights below the plane 
for each radial cut. We note that stars with a distance of about 6 kpc 
on the far side of the bulge at b — —10° are located at z— heights 
of 2.4 kpc below the plane. At b = -10° and \R G \ > 7 kpc the 
stars are on the far side of the bulge at mean z— heights > 2.6 kpc. 
Therefore it is not surprising that the component C (with its thick- 
disk-like metallicity range) should dominate on the far side of the 
bulge. 



4.2 The metallicity gradient with distance \ z\ from the plane 

Along the minor axis, the mean [Fe/H] for some of the individual 
components of the MDF shows a weak decrease with latitude. From 
b — — 5° to b = —10°, the total decrease in mean abundance 
of components A and B is -0.06 ± 0.04 dex and -0.06 ± 0.03 
dex respectively. The mean abundance of component C shows no 
significant decrease with latitude on the minor axis, although the 
longitude-integrated MDF does show a small decrease in the mean 
of component C with latitude (see Table [3} of —0.03 ± 0.02 dex. 
The numbers of stars with [Fe/H]< —1 is small (10, 29 and 55 
stars at b — —5°, —7.5° and —10° respectively along the minor 
axis) and we cannot evaluate the abundance trends with height for 
the metal-poor components D and E. 

A small decrease in the mean abundances of components A 
and B for the integrated MDFs (Z = ±15°) is also measured for 
these components. This small change in the mean abundances with 
latitude, and the large change of the contribution fractions of the 
components, indicates that the overall observed abundance gradi- 
ent measured by Zoccali et al. (2008) of about —0.6 dex kpc -1 
from b — —4° to b — —12° comes mainly from the changing pro- 
portions of the MDF components with latitude, as shown in Figure 
II II rather than from gradients in the individual components. 

From b = —5° to b = —10°, the gradients for A and B with 
height below the plane calculated from the minor axis MDFs are 
-0.08 ± 0.04 dex kpc -1 . No gradient in C is seen along the minor 
axis but the integrated MDFs show a small gradient of — 0.04±0.02 
dex kpc - . 



4.3 The radial metallicity gradients 

The mean abundances for components A, B and C as a function of 
Rg are given in Table [5] with their corresponding errors. A weak 
radial gradient (~ —0.02 dex kpc -1 ) is evident in the mean abun- 
dance of the most metal-rich component A. We note that this gra- 
dient may come partly from a radial gradient and partly from a 
z— gradient in these zones of constant latitude. This weak abun- 
dance gradient is shallo wer than the gradient seen in the cepheids 



Luck & Lambert 



Yong. Carney. & Friell 



20 1 1 ) and young (age ^ 1 Gyr) open clusters 
20121) . However, given an additional neg- 



ative gradient with height from the plane, this measured gradient 
will be smaller than the true value. This is because, with increasing 
distance along the line of sight, the stars sampled are further from 
the plane and will be relatively more metal poor. 

At Rg > 8 kpc at low latitudes and Rg > 4.5 kpc at high 
latitudes, the most metal rich fraction of stars is no longer present 
and the MDF can be described by two main components: the more 
metal rich component is component B. This component represents 
the most metal rich stars that we sample at large heights from the 
plane; we note that at b — —5° the |z|-height at 6 kpc on the far side 
of the bulge is 1.2 kpc. From the inner region, component B shows 
almost no abundance gradient with | Rg | . The abundance of com- 
ponent B is higher on the far side of the bulge where components A 
and B have merged (b = — 5° ) or component A is no longer present 
(b = -10°). 

The mean abundance of component C remains relatively flat 
at both latitudes out to about \Rg\ = 4.5 kpc. Outside this radius, 
the errors are too large to evaluate any gradient. At _Rg > 8kpc, the 
I z\— height above the plane is about 3 kpc for b = —10° and there 
may be increasing contamination of component C by the metal- 
poor stars which belong to the overlapping lower metallicity com- 
ponent D. 

We have used the above gradient in component A to argue in 
the previous section that components A and B merge at a radius 
\Rg\ « 4.5 kpc. Within R G < 4.5 kpc at both latitudes, the MDF 
shows three significant peaks; two at [Fe/H] > -0.5 and one at -1.0 
> [Fe/H] > -0.5. Outside of Rg > 4.5 kpc, there is only one clear 
peak in the MDF present at [Fe/H] > -0.5 and a second significant 
peak at -1.0 > [Fe/H] > -0.5. At all radii there is the small con- 
tribution of stars with [Fe/H] < -1.0 represented by components D 
and E. 



5 ALPHA-ENHANCEMENT 

The a— enhancement reflects the rate of chemical evolution. Stars 
are a— enhanced if the chemical enrichment in their formation en- 
vironment has been sufficiently rapid that there was not time for 
the Fe-enriching SN la to contribute fully to their chemical evolu- 
tion. In the following, [a/Fe] denotes the mean of [Si/Fe], [Mg/Fe] 
and [Ti/Fe] (we used 5 Si lines, 3 Mg lines and 9 Ti lines); we 
were unable to make an accurate estimate of [Ca/Fe]. Figure IT61 
compares the radial distribution of a— enhancement for stars in 0.5 
dex bins of [Fe/H] from -2.0 to 0.6, extending in radius from in- 
side the bulge to the outer disk and halo at Rg = 30 kpc. For this 
figure, we have not associated the [Fe/H] intervals with the MDF 
components, because of the abundance gradient seen in the most 
metal-rich component A. 
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Table 5. Mean metallicity of Gaussian components A, B and C at b = —5° and b = —10°, for \l\ < 15°, as a function of Galactocentric radius, \Rq |. The 
range of \z\ -height from the plane for each range is indicated. 



Distance cut 


Number of 




Mean [re/HJ 








stars 










\Ro\ (kpc) 


b=-5° 


\z\ kpc 


A 


B 


C 


<0.75 


300 


0.6-0.8 


0.16±0.02 


-0.24 ±0.02 


-0.67±0.04 


<1.5 


1100 


0.6-0.8 


0.15±0.01 


-0.22±0.01 


-0.7 ±0.03 


1.5-3 


2300 


0.1-0.4 


0.13±0.01 


-0.27±0.01 


-0.69±0.03 


3^.5 


2000 


0.1-0.3 


0.10±0.01 


-0.28±0.01 


-0.67±0.02 


4.5-6 (near) 


600 


0.2-0.3 


0.1 1 ±0.03 


-0.26 ±0.04 


-0.75±0.07 


5-7 (near) 


200 


0.1-0.3 


0.11 ±0.03 


-0.28 ±0.04 


-0.67±0.07 


4.5-6 (far) 


120 


1.1-1.2 


0.12 ±0.08 


-0.25 ±0.06 


-0.70±0.07 


5-7 (far) 


100 


1.1-1.3 


0.10 ±0.08 


-0.25 ±0.06 


-0.68±0.07 


8-10 (far) 


70 


1.4-1.6 


- 


-0.13 ±0.04 


-0.75±0.08 


10-15 (far) 


130 


1.6-2.0 


- 


-0.26 ±0.08 


-0.77 ±0.08 


15-40 (far) 


100 


2.0-4.2 


- 


-0.21 ±0.08 


-0.95 ±0.07 


\Rq\ (kpc) 


b=-10° 


\z\ kpc 


A 


B 


C 


<0.75 


180 


1.3-1.5 


0.13±0.06 


-0.27±0.03 


-0.67±0.04 


<1.5 


850 


1.1-1.7 


0.10±0.03 


-0.28±0.02 


-0.68±0.01 


1.5-3 


2450 


0.9-1.9 


0.09±0.03 


-0.29±0.02 


-0.65±0.01 


3^.5 


1700 


0.9-2.2 


0.08±0.03 


-0.27±0.02 


-0.67±0.01 


4.5-6 (near) 


250 


0.3-0.6 


0.10 ±0.06 


-0.22 ±0.05 


-0.71±0.05 


5-7 (near) 


150 


0.2-0.5 


0.04 ±0.05 


-0.25 ±0.05 


-0.76±0.04 


4.5-6 (far) 


200 


2.2-2.5 




-0.35 ±0.05 


-0.78±0.06 


5-7 (far) 


120 


2.3-2.6 




-0.25 ±0.08 


-0.77±0.06 


8-10 (far) 


40 


2.8-3.1 




-0.27±0.08 


-0.79±0.08 


10-15 (far) 


150 


3.1^.0 




-0.13±0.08 


-0.91±0.06 


15-40 (far) 


170 


4.0-8.4 




-0.20±0.08 


-0.70±0.06 



5.1 Alpha-Enhancement trends with Galactic radius 

Stars with [Fe/H] > —0.5, and particularly the most metal rich 
stars, show a decrease in [a/Fe] from slightly enhanced values in 
the inner Galaxy. Stars with [Fe/H] > show [a/Fe] declining to 
negative values beyond \Rg\ > 20 kpc (Figure |T6j. This decrease 
in the [a/Fe] ratio with radius indicates that the chemical evolution 
in the inner Galaxy proceeded more rapidly than in the disk at larger 
radii where the [a/Fe] ratio is near solar for the most metal-rich 
stars. Note that stars at these larger values of \Rg\ > 7 kpc are 
mostly on the far side of the bulge, at z— heights > 1.4 kpc. 

For stars with —2.0 < [Fe/H] < —0.5, the a— enhancement 
is fairly flat with galactic radii out to 30 kpc. The most metal poor 
fraction of stars shows a slight increase in enhancement, although 
these stars show a larger scatter. 

Figure [T7] shows the radial distribution of mean [a/Fe] for the 
most metal-poor stars with [Fe/H] < -2.0. (The most metal poor 
star in the sample is at [Fe/H] = —2.9; however it is difficult to 
measure the gravity from the weak Ca-triplet lines and we have no 
distance estimate for this star). The distribution of [a/Fe] for these 
metal-poor stars show a steep decrease from the inner region. In this 
subsample of metal-poor stars with |i?o| < 8 kpc, 75% lie within 
\R G \< 4.0 kpc. 

5.2 Alpha-enhancement and components 

We can analyse the a— enhancement in the context of the com- 
ponents of the MDF Figure Q~8] shows contour plots of the stellar 
density for stars with Rq < 3.5 kpc in the [Fe/H] - [a/Fe] plane, 
within longitude ±15°, with the generalised MDFs for the compo- 
nents ABC superimposed. At 6 = —5° there are two clear peaks 
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Figure 16. Alpha-enhancement of stars in 0.5 dex intervals of [Fe/H] 
against Galactic radius | Rq \ . The size of the radial bins increase with in- 
creasing distance. 



in the density distribution, corresponding to components A and B. 
At b — —10°, component A has attenuated relative to component 
B and the contour plots show only a single density maximum asso- 
ciated with component B. The most metal-rich component A has a 
mean [a/Fe] of about 0.1 at b = -10° and 0.16 at b = -5° (in 
agreement with the weakly enhanced [a/Fe] at solar [Fe/H] found 
by high resolution stud ies of bulge stars in Baade's window by 
lAlves-Brito et aljfeoiol) '). All of the lower-metallicity components 
in the bulge are more a— enhanced. The a— enhancement increases 
as expected with decreasing [Fe/H], from about 0.1 at solar metal- 
licity to about 0.6 at [Fe/H] = —2.0. The dispersion in [a/Fe] in- 
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Figure 17. Alpha-enhancement of our 27 most metal-poor stars with [Fe/H] 
<-2.0 against | Rq \ . For | Rq | > 8 kpc, we have only three metal-poor 
stars (at \Rq\ ~ 9kpc, llkpc and 12 kpc). 




-1.5 -1.0 -0.5 0.0 0.5 
[Fe/H] 



creases as [Fe/H] decreases. At latitude —10°, the contribution of 
the more metal-rich component of stars is attenuated, and we see a 
lower dispersion in the [a/Fe] for [Fe/H] > relative to the lower 
latitude strip. The \ 2 error on the [a/Fe] measurement is a = 0.1 
and this is consistent with the spread seen in Figured! 

Figure [T9] compares the [Fe/H] - [a/Fe] relation for stars in 
three radial zones: inside the bulge at \Rg\ < 1-2 kpc and \Rg\ < 
3.5 kpc, and outside the bulge at 5 kpc < \Rg\ < 10 kpc. For the 
innermost interval inside the bulge region, the a— enhancement in- 
creases linearly with decreasing [Fe/H], from about 0.1 for stars 
above solar metallicity to about 0.6 at [Fe/H] = —2.0. The disper- 
sion in [a/Fe] increases more slowly, from about 0.1 at the most 
metal-rich end to about 0.15 at [Fe/H] = —2.0. This dispersion at 
the more metal-poor end is higher for stars within \Rq | < 3.5 kpc, 
at a ~ 0.35; other dispersion values in this radial interval follow 
those of the inner region Rq < 1 kpc. 

Comparing the a— enrichment of stars in the disk and the 
bulge, we see that the stars with 5 kpc < | Rq | < 10 kpc and [Fe/H] 
> —0.5 have similar [a/Fe] values to those inside the bulge. The 
[a/Fe] values are about 15% lower for stars with [Fe/H] < —1.0 
and about 10% higher for stars with [Fe/H] > -1.0 relative to those 
with |_Rg| < 3.5 kpc. The dispersion is higher than for \Rg\ < 3.5 
kpc by about 10% at the more metal-rich end and 30% at the more 
metal-poor end. 



6 INTERPRETING THE COMPONENTS 

In this section, we discuss two significant features of the MDF com- 
ponents: (1) their individual vertical abundance gradients and the 
overall vertical abundance gradients associated with the changing 
fraction of the major components with height above the plane, and 
(2) the contribution of the components to the split red clump of the 
bulge. 



6.1 The vertical metallicity gradient of the bulge 

In the bulge region, components A, B and C individually show 
weak vertical abundance gradients (section 4.2) of —0.08 ± 0.04 
to -0.04 ± 0.02 dex kpc -1 over the latitude interval b = -5° 
to b — —10°. Although our overall vertical metall icity gradient 
(—0.45 dex kpc -1 ) is comparable to that found by IZoccali et al.l 
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Figure 18. Alpha-enhancement at b = —5° (top) and b = —10° (bot- 
tom) in the bulge region RgI < 3.5 kpc. The MDF components at each 
latitude are superimposed for comparison. The grayscale colour legends at 
right show the number of stars. 



d2008l) , we interpreted this overall gradient as coming from chang- 
ing fractions of components A-D with latitude via the action of the 
instability process on the pre-existing early disk. 

It is not clear how to interpret these small vertical abundance 
gradients in the individual components. Are they the remnants of 
abundance gradients that were already in place in the thin and thick 
disks before the instability event? The small gradients which we 
observe compar e well with the fairly flat ve rtical abundance distri- 
bution found by Bek ki & Tsuiimotol l l201ll) . from their simulation 
of an initial disk system which buckles into a boxy/peanut bulge/bar 
via the disk/bar instability. They model both the instability of an ini- 
tial single thin disk and also a dual (thin + thick) disk population. 
These disks initially have vertical abundance gradients. Due to the 
vertical mixing of the stellar populations by the bar, they find that 
the bulge is unlikely to have a steep vertical metallicity gradient. 
Their model does not include chemical evolution so is appropri- 
ate for a system in which the chemical evolution took place in the 
early disk, before the instability events and the formation of the 
bar/bulge. 
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Figure 19. Mean [a/Fe] and dispersion in [a/Fe] versus [Fe/H] for three 
intervals of Rq: in (a) Rg < 1.2 kpc, and (b) < 3.5 kpc, the stars lie 
within the bulge region, and in (c) 5 kpc < Rq < 10 kpc they lie outside 
the bulge region. 



6.2 The split red clump of the Galactic bulge 

There is a striking relationship between the spatial distribution of 
stars in the bulge region and [Fe/H] which was seen in our analysis 
of the split clump (Ness et al. 2012a). The split clump is a conse- 
quence of the peanut structure of the bulge, so the presence of a 
split clump indicates membership of the peanut bulge. 

Figure l20l shows the magnitude distribution of the red clump 
stars along the minor axis for our three latitudes, selected with Ko = 
12.38 to 13.48 and log g = 1.9 to 3.1. The clump stars are split into 
three [Fe/H] bins; (i) [Fe/H] > 0, (ii) > [Fe/H] > -0.5 and (iii) 



[Fe/H] < -0.5. In Ness et al. (2012a) we compared the spatial dis- 
tribution and kinematics of the clump stars with predictions from 
an evolutionary N-body model of a bulge that grew from a disk 
via bar-related instabilities. The density distribution of the peanut- 
shaped model is depressed near its minor axis, thus producing a bi- 
modal distribution of stars along the line of sight through the bulge 
near its minor axis, very much as seen in our observations for stars 
[Fe/H] > -0.5. We found that only stars with [Fe/H] > -0.5 show 
the split structure at latitudes b < —5°, and concluded that the stars 
of the boxy /peanut bulge have [Fe/H] > -0.5. Stars with [Fe/H] < 
-0.5 do not show the split clump at any latitude. This indicates that 
the initial disk from which the bulge formed had very few stars with 
[Fe/H] < -0.5, or that stars with [Fe/H] < -0.5 were not in a region 
of phase space which is mapped into the boxy structure by the disk 
instability. 

From Figure [20] stars with [Fe/H] > -0.5 show the split at 
b < —5°. Stars in our component A ([Fe/H] > 0) show a deeper 
split than stars in our component B (0 > [Fe/H] > —0.5) at the 
intermediate latitude of b = —7.5°. At b = —10° the split is still 
seen for components A and B but is affected by the asymmetry in 
the magnitude distribution for the more metal-rich stars: these stars 
are located preferentially closer to the plane and are therefore found 
more prominently in the nearer and brighter peak. The proportion 
of stars with [Fe/H] > decreases from about 40% at b = —5° to 
about 8% at b — —10°. In contrast, for the component B stars with 
—0.5 < [Fe/H] < 0, the proportion stays approximately constant 
from b = -5° to b = -10°. 

The split distribution of the component A and B stars, and the 
absence of the split for the more metal-poor component C stars, 
leads us to infer that the peanut bulge structure is defined primarily 
by the stars of component A ([Fe/H] > 0) and component B (0 > 
[Fe/H] > -0.5). 



7 DISCUSSION 

In a survey of the inner Galaxy, we can expect to find stars from 
the inner regions of all of the major components of the Galaxy. 
We have identified five metallicity components in total in the 
stellar MDF of the inner Galaxy, and interpreted the MDF as a 
composite population made up of these components A-E in order 
of decreasing [Fe/H]. The two more metal-rich components (A 
and B) define the boxy/peanut structure of the bulge. We make the 
following interpretation of the metallicity components within the 
disk instability scenario for bulge formation: 

• A: The metal-rich "boxy/peanut-bulge" with a mean [Fe/H] ~ 
+0.15 that is concentrated towards the plane. In the inner regions, 
component A has the chemical properties of a younger chemically 
evolved thin disk (higher [Fe/H] compared to the thin disk in the 
solar neighbourhood and slightly higher [a/Fe]). At larger Galactic 
radii, component A appears to be somewhat more metal-poor (see 
Table 5) and we associate it with the thin disk at larger Galactic 
radii. Its MD F at larger radii is similar to that of the thin disk near 
the Sun (e.g. lHavwoodll2008h . Although chemically similar to the 
thin disk, it replicates the line-of-sight kinematics and spatial prop- 
erties of the thicker boxy /peanut-bulge com ponent B in the inner 
regions of the Galaxy (see lNess et alj2012cl in preparation), but at 
lower heights above the plane (see Section 4). We note again that 
stars in the metallicity range of component A are prominently in- 
volved in the bimodal peanut structure of the bulge seen at higher 
Galactic latitudes. 
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Figure 20. The magnitude distributions for the red clump stars selected with Kq = 12.38 to 13.48 in three [Fe/H] intervals of (i) [Fe/H] < -1.0, (ii) -0.5 < 
[Fe/H] < , (iii) [Fe/H] > 0. Gaussian smoothing kernel a = 0.15 mag was applied. The y-axis is the normalised number of stars. 



• B: The vertically thicker "boxy /peanut-bulge" contributes 
a similar fraction of the bulge stars over the latitude range from 
b = -5° to & = -10°. Its mean [Fe/H] value is about -0.25 
and it is the dominant component in most of our bulge fields inside 
Rg < 4.5 kpc (see Figures 10 and 11). In the instability scenario, 
it would represent the early thin disk when there were few stars 
with [Fe/H] < — 0.5 dex at the onset of the bar buckling instability. 
Again, stars in the metallicity range of component B show the split 
structure of the bulge at higher latitudes. At larger radii, component 
B is associated with the metal poor stars of the early thin disk that 
we sample at large heights from the plane. 

• C: We interpret component C as the inner thick disk (mean 
[Fe/H] about -0.70) which is at most weakly involved in the 
boxy/peanut structure of the bulge. Component D is tentatively 
ascribed to the metal-weak thick disk (mean [Fe/H] about —1.2). 
These components were likely present in the inner regions of the 
Galaxy (] Rq\ < 3.5 kpc) before the instability event, and may have 
been puffed up by this event. We note again that components C and 
D do not appear to be part of the split structure of the peanut bulge 
(Ness et al., 2012a). At larger j Rg\, these old components are still 
present in their relatively undisturbed form. 

• E: We interpret this metal-poor component as the stars of the 
inner Galactic halo. 

Components A and B are clearly part of the boxy/peanut struc- 
ture of the Galactic bulge, as show n by the split red clum p discussed 
in Ness et al, 2012a. Models by lAthanass oula ( 2005]) show dif- 
ferent levels of boxy/peanut behaviour. For the boxy systems, the 
isophotes are roughly parallel to the equatorial plane. The isophotes 
of the peanut systems show a clear dip near Rg = 0. If this dip 
is relatively deep, the systems give an impression of an X-shaped 
structure most clear at larger heights above the Galactic plane, par - 
ticularly after unsharp masking. We have seen in lNessetal.lf2012j) 
that component A shows a more prominent central minimum than 
component B, but it is not clear whether A is intrinsically more X- 
shaped or if its deeper minimum reflects a similarly shaped struc- 
ture but with a smaller scale height than the thicker boxy/peanut- 
bulge represented by component B. 

7.1 The origin of the components 

In the instability scenario, the boxy/peanut-bulge formed out of the 
early thin disk: after about 1 Gyr, the instability of the disk forms 
a bar which then buckles vertically into the observed boxy/peanut 
morphology. This scenario gives a guide to interpreting the compo- 
nents of the bulge. We know from the age-metallicity relation in the 



solar neighbourhood (e.g. Haywood 2008) that the chemical evo- 
lution of the solar neighbourhood proceeded rapidly: 10 Gyr ago, 
the stars of the thin disk in the solar neighbourhood already had a 
mean [Fe/H] near solar, with a spread from about +0.5 to —0.5. 
In the context of inside-out formation of the disk, it is likely that 
the chemical evolution of the inner Galaxy was even more rapid. 
By the time that the bar-forming instability occurred, the chemical 
evolution of the inner disk was probably already well advanced. 

The instability process would affect all of the inner disk stars, 
and generates a mapping of the stars of the early disk into the 
boxy/peanut structure. In this way, the bulge preserves a dynam- 
ical imprint of the chemical distribution of the thin disk at the time 
that the buckling occurred; the bulge is a chemical snapshot of the 
MDF of the early thin disk captured in the extended bar. We have 
seen that two components, A and B, are involved in the boxy/peanut 
structure of the bulge. Figure QT] shows that the more metal-rich 
component A (mean [Fe/H] = +0.15) is more concentrated to the 
plane than the thicker dominant bulge component B (mean [Fe/H] 
= —0.25). This more metal-rich component A is chemically con- 
sistent with what we would expect for a younger, colder and more 
metal-rich thin disk in the inner region at the time of the instability 
event. 

Examination of the N-body models shows that the mapping 
from disk to bulge depends on the location and kinematics of the 
stars at the time of the instability. Kinematically colder stars in the 
disk can suffer substantial radial and vertical migration, and can 
therefore be strongly involved in the boxy/peanut structure. We will 
discuss this mapping of the disk into the bulge in more detail in a 
later paper. At this stage, we suggest that component A may simply 
represent the responsive cold and metal-rich stars of the thin disk in 
the inner region which were entrained in the potential of the major 
boxy/peanut component B. The kinematical selectivity of the map- 
ping is likely to explain the presence of the two distinct components 
A and B in the bulge at the higher latitudes | b \ > 5° . 

The thick boxy/peanut bulge population B can be interpreted 
as the product of the bulk of the thin disk (already a few Gyr old 
at the time of the instability event) which form s a bar and puffs 
up into a boxy/peanut-bulge ( Athanassoula 2008). The mean abun- 
dance and abundance range of component B is similar to that of the 
intermediate-age disk stars (3 to 8 Gyr old) in the solar neighbour- 
hood at the present time (e.g. Haywood 2008). The bulge structure 
is visible in NIR photometry and star counts out to a radius Rg of 
about 3.5 kpc. Further from the Galactic center, stars of the old thin 
disk (component B) are still present, but in their relatively undis- 
turbed original thin disk structure. 
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Component C has a mean abundance like that of the thick disk 
in the solar neighbourhood. If it is representative of the thick disk 
in the inner Galaxy before the instability event, it may have been 
correspondingly hotter and less dynamically responsive at the time 
of the instability. That may be the reason why component C does 
not appear to be involved in the boxy/peanut structure of the bulge. 
Similarly, component D has a mean [Fe/H] of ab out —1.20, like 
that of the metal- weak th ick disk near the Sun (e.g. lChiba & Beersl 
l200CtlCarollo et alfcOlOh . 

Component E with mean [Fe/H] = —1.7 appears in small 
fractions at about the same [Fe/H] in most fields, and shows no ev- 
idence for an abundance gradient with latitude. It is slowly rotating 
(see Ness et al., 2012c) and is probably an entirely separate popula- 
tion, not kinematically or chemically associated with the disk. We 
interpret the metal-poor component E stars as inner halo stars on 
orbits of high eccentricity passing through the bulge region. They 
may correspond to the high eccentri city stars at an abunda nce of 
about —1.7 i dentified near the Sun bv[chiba & Beersl d2000|). Sim- 
ulations (e g.lDiemand. Madau. & Moorell2005MBrook et alj|2007l ; 
lTumlinson]|2010l) show that the oldest stars in the Galaxy, formed 
in fragments from high-a density fluctuations, should now inhabit 
the inner Galaxy. In this context, the stars of component E in the 
inner Galaxy would be worth further spectroscopic study at higher 
resolution. 

We now return to the metal-rich component A, which is 
more prominent closer to the Galactic plane but at the same time 
shows most strongly the split red clump which we associate with 
the boxy/peanut structure of the bulge. The MDF of the solar 
neighbourhood may again provide a useful reference. In the solar 
neighbourhood, the youngest disk stars (ages < 2 Gyr) are rela- 
tively metal-rich, have a relatively narrow range of abundance, and 
are kinematically colder than the intermediate-age disk stars. The 
larger width of the MDF for the intermediate-age stars in the solar 
neighborhood is currently attributed to the secular effects of radial 
migration, which radially smooths the Galactic abundance gradi- 
ent and so increases the width of the local MDF. In the same way, 
the more metal-rich stars of the early inner thin disk were prob- 
ably kinematically colder and dynamically more responsive, thus 
preferentially mapped into boxy/bulge-supporting orbits. 

Alternatively it is possible that component A comes, via 
a second later and less vertically extended bar-forming insta- 
bility event, from the more metal-rich thin disk which had 
continued to undergo star formation and chemical evolution 
since the first instability episode. Simulat ions have shown that 
a second buckling episode can occu r dAthanassould l2005bl : 
iMartinez-Valpuesta. Shlosman. & HelledbOOrj) . well after the first 
one, and that it leads to a vertically thicker structure than given by 
the first only buckling event. If we make the plausible assumption 
that component B forms earlier from the older metal-poor disk, and 
component A forms later from the metal-rich disk, then component 
B may well have had sufficient time to undergo a second buckling, 
while component A has not. This offers a possible explanation for 
why component B is more vertically extended and more metal poor 
than A. A secondary instability could also be triggered by a sig- 
nificant accretion event. Models with two disks of different scale 
heights may also offer more insight. The current N-body models 
we are working with do not show a secondary instability event and 
have a single disk only; we therefore cannot test these effects. 

Component A attenuates in density by a factor of 5 between 
b — —5° and b = —10°. It is possible that the original bar 
formed in the early disk may have gradually entrained the chem- 
ically evolving thin disk during and after the instability process. 



In this way, the more metal-rich stars could have experienced only 
the later fraction of the bar-buckling process, and remained more 
closely bound to the Galactic plane. 

In these alternative scenarios, we are associating component 
A with the inner thin disk which has continued to undergo star for- 
mation and chemical evolution. We can trace this component out to 
large Galactic radii, where it becomes gradually more metal-poor 
(Table 5), reflecting the abundance gradient of the Galactic disk. 
In the central region however this component behaves in a clearly 
bar-like way, reflected in the split clump (see Section 6.2) and the 
kinematics as a function of longitude (Ness et al., 2012c). 

7.2 The alpha-enhancement of the components 

Fi sure [T8l shows that the [a/Fe]-[Fe/H] relation for stars in the in- 
ner Galaxy follows the pattern in the solar neighbourhood, except 
that the most metal-rich stars in the inner regions do show a weak 
a-enhancement. Alpha-enhancement is a measure of the rate of 
chemical evolution. The more metal-poor boxy/peanut component 
B, which we associate with the older pre-bulge disk, is a-enhanced 
above the level of the more metal-rich bulge component A which 
we associate with the younger pre-bulge thin disk. This is consis- 
tent with our argument that the bulge stars of component B were 
formed earlier as part of the older disk. Their a-enrichment would 
be present not only in the bulge population but also for the other 
stars of the early disk which are still now in the disk. 

The more metal-poor components in our survey, with [Fe/H] 
< —0.5, include the thick disk and metal-weak thick disk. Al- 
though there are several possibilities for the formation of the thick 
disk, the stars of components C and D may also have originated 
from the even earlier thin disk, heated kinematic ally by early mi - 
nor merger activity before the instability event (e,g. lFreemanl2008h . 
They would be the oldest stars of the original disk, and their level 
of a— enhancement indicates that the duration of their star forma- 
tion was short. Some of these more metal-poor stars may also have 
been puffed up into the bar/bulge by the instability event itself. The 
a-enhancement for all stars with [Fe/H] < -0.5 is fairly flat out to 
Rg = 30 kpc. As their a— enhancement is directly related to the 
duration of their star formation event, this suggests that the more 
metal-poor stars which contribute to the bulge/bar, the thick disk 
and metal-weak thick disk may have had a common origin in the 
oldest stars of the initial Galactic thin disk. In this picture, the dif- 
ferent MDF components are fossil remnants from the dynamical 
redistribution of an initial disk system of the Milky Way. The abun- 
dances and a-enhancement are proxies for the time at which the 
stars of each component formed and the rate of their formation. 
Chemodynamical models may be able to link the abundances and 
timescales for this formation and structural dispersion. 

Figure [T8] shows the clearly bimodal distribution of abun- 
dances between the two boxy/peanut components A and B in our 
b = — 5° fields. Component A has a mean [a/Fe] value of 0. 16, and 
the old disk component B has a mean [a/Fe] of 0.23. Several au- 
thors ha ve argued that the bulge is an a-enhanced metal-rich pop- 
ulation feoccali et alj|2008t iFulbrigtit. McWilliam. & Richll2007l ; 



iLecureur et al .120071 : iMelendez et alj2008h . We see that the metal- 
rich component A is only weakly a-enhanced, and the thicker 
a— enhanced bulge component B is not so metal-rich. 

From Figure [19] we can compare the a-enhancement of stars 
with similar [Fe/H] abundances in the bulge and in the disk at larger 
Rg - For stars with [Fe/H] > — 1, which we would identify as origi- 
nating from the disk, the mean [a/Fe] values are very similar in the 
bulge (Figures[l9](a) and (b) ) and the disk (Figure [19] (c) ). For the 
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more metal-poor stars with [Fe/H] < — 1.0 the alpha-enhancement 
at large |_Rg| (Figure [T9](c)) is similar to stars in the inner region 
\Rg\ < 3.5 kpc. However the dispersion in the alpha enhancement 
is significantly lower in the inner bulge region (Figure |19b com- 
pared to at larger radii. This may reflect that the metal-poor popula- 
tion in the centre is more homogenous than in the outer regions. For 
the very metal-poor stars with [Fe/H] < — 2, the a-enhancement is 
significantly higher in the inner bulge region than in the disk at 
larger radii (see Figure [771). indicating a more rapid star formation 
history for these very metal-poor stars in the inner regions of the 
Galaxy. 



8 COMPARISON TO OTHER STUDIES 

We have found an overall metallicity gr adient in the whole MDF 
across our latitudes, comparable to that o f lZoccali etal] <f2008). Our 
MDF at b — —10° and our a-enhancement res ults at our higher 
latitud e, also compare well with the results of lUttenthaler et all 
d2012h for their high resolut ion study at (I, b) = (0°, -10°). Note 



Uttenthaler et alj J2012h also confirm the presence of the split 



that 

clump only for the more metal rich stars. Other spectroscopic stud- 
ies of bulge stars have found evidence for multiple components 
in the MDF, and we now compare our results with theirs. Ta- 
ble [6] summarises the mean abundances for the components found 
in the ARGOS survey compared with tho se f rom the work b y 
iBabusiaux etafl j201of ). lBensbv et al.1 d20ld) and lHill et all d201ll) . 
Each of these studies sees some but not all of the components found 
in the larger ARGOS sample. 

8.1 Vertex deviations and the bulge 

In the study of about 400 bulge K-giants and clump giants in 
Baade's win dow at b = —4 ° and a field near the minor axis at 
b = -6° bv lBabusiaux et alj l l20icl) . our components A, B and D 
can be seen. They argue that both of the main scenarios of bulge 
formation have contributed to the bulge. Using proper motion and 
radial velocity data, they measure the vertex deviation of the ve- 
locity ellipsoid for their identified components. Their more metal- 
rich component, with [Fe/H] near 0.14 (our A), shows a signif- 
icant vertex deviation and they conclude that this component is 
the bar/bulge that formed via the instability scenario. Their more 
metal-poor components (our B and D) show a smaller vertex de- 
viation, and they associate these components with a spheroid and 
thic k disk. 

IBabusiaux et alj ( feoioh also measured the radial velocity dis- 
persion as a function of [Fe/H] for their fields at b = —4°, —6° 
and —12° near the minor axis of the bulge. They find that the dis- 
persion increases with [Fe/H] in their b — — 4° field (the more 
metal-rich stars are kinematically hotter) but decreases as a func- 
tion of [Fe/H] in their higher latitude fields. Similarly, we find that 
at b — —5°,— 7.5° and —10°, the velocity dispersion decreases 
with [Fe/H] along the minor axis (and indeed at all longitudes, see 
iNess et al]|2012cl in preparation). This change with latitude of the 
sense of the abundance - (velocity dispersion) relationship indicates 
a structural change of the bulge populations at around b — —4°. 

We now discuss the vertex deviation measurements of the 
bulge near the minor axis and compare these results to our N- 
body model. We do not have proper motions for our stars so 
cannot independently measure the vertex deviations. The vertex 
dev iation of the stars in Baade 's win dow has been measure d 
by ISoto. Rich. & Kuiikenl d2007l) and IBabusiaux etail d20ld) . 



ISoto. Rich. & Kuiikerj J2007h conclude from their measurements 
that stars with [Fe/H] > —0.5 are supporting a bar. They find no 
significant vertex deviation for stars with [Fe/H] < —0.5. This 
makes sense in the context of our finding of the metallicity depen- 
dence ofjte_splk_chrm2_and assignment of components to popula- 
tions. [Soto[jijch^&iKui]kej3 d2007l) find a trend of decreasing ve- 
locity d ispersion with in c reasin g metallicity as do lBabusiaux et all 
d2010h . I Babusiaux ct al. I d2010h also find the largest vertex devia- 
tion for the more metal-rich stars which we would associate with 
component A. Their vertex deviation in the Ir plane (denoted l v ) for 
stars with [Fe/H] < —0.14 is —13° ± 9. For the more metal-rich 
stars with [Fe/H] > —0.1 the vertex deviation is about —40°. They 
argue that the most metal-rich stars represent the bar/bulge in the 
central region, and the stars with lower metallicities are associated 
with an old spheroid. 

We consider first the vertex deviation of stars belonging to 
our boxy/peanut-bulge component B, with > [Fe/H ] > — 0.5. 
The l v vertex deviation measured by IBabusiaux et al.1 d20ld) for 
stars in this abundance range is about —30°. This agrees well with 
the estimate from our bar/bulge model, for which we find a vertex 
deviation l v = —26° ±2° for all stars within a galactocentric radius 
of 0.5 kpc. 

Now we address the vertex deviation of stars belonging to our 
boxy/peanut-bulge component A, with [Fe/H] between and 0.5. 
For stars in this abundance range, IBabusiaux et al.1 d2010h measure 
a vertex deviation in Baade's window of l v — —43° ± 5° . We 
examined 6 models with different shapes and scale heights of the 
central bulge, scaled to match the kinematics of the Milky Way, 
and found the vertex deviation near the minor axis to be fairly in- 
sensitive to the scale heights of the bulge/bar and also to the height 
above the plane. To reproduce this higher l v in our models, we 
needed to rotate the bulge/bar in the model from the previously 
adopted 20° to 45° with respect to the line of sight. The vertex de- 
viation in the model for a bar angle of 45 ° is l v = —45° ± 2° 
similar to that derived by IBabusiaux et al.1 J2010l). This is an in- 
teresting result in light of the work oflCabrera-Lavers et al.1 ( l2008h 
and others (e.g. IWeinberg|[T992l:lHammerslev et alj|l994^ who ar- 
gue for the existenc e of a long bar in the Milky Way from analy- 
sis of 2M ASS data dSkrutskie et al]|2006l) . Using red clump giants 
as tracers. Icabrera-Lavers et a" 1 20081) present evidence for a long 
bar with a semi-length of 4 kpc, ending around I — 28°. This long 
bar inclined at 45° coexists with the thicker bulge/bar inclined at 
20° and ~ 3.5 kpc in length dGerharj|2002h (our component B) . 
It has been previously suggested (e.g. ICabrera-Lavers et al. I l2008h 
that the inner Galaxy is dominated at higher latitudes by the bulge 
and at lower latitudes by the long bar. 

We are not suggesting that the vertex deviation in our model 
is strong evidence for the coexistence of a boxy/peanut-bulge made 
of stars with > [Fe/H] > —0.5 and a long flatter bar that con- 
tains the metal-rich stars with [Fe/H] > in the bulge. It is however 
an interpretation, in light of the dependence of vertex deviation of 
the model on the bar orientation. Changing the angle however may 
be one of several ways to affect the vertex deviation. Given the 
magnitude distribution of the split clump at the higher latitudes as 
a function of [Fe/H] and the changing relationship between dis- 
persion and [Fe/H] at b = —5°, it seems clear that there are two 
components to the bulge population. It is possible that these results 
could be explained by a secondary structu r e whic h is dominant at 
low latitudes as per ICabrera-Lavers et ai] d2008l) . The metal-rich 
stars may be associated with a thin disk population associated with 
a bar orientated at 45° to the Sun-centre line. 
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8.2 Microlensed dwarfs in the bulge 

A potential problem with our interpretation o f component B as the 
primary component of the bulge comes from lBensbv et al.l feOllh 
who compiled abundances for 26 microlensed dwarfs in the direc- 
tion of the bulge. The [Fe/H] values for their stars extend from 
—0.72 to +0.54. The MDF of their microlensed dwarfs shows two 
well-separated components corresponding to our components A 
and C (which we argue are the thin boxy/peanut-bulge and the old 
thick disk). As for our sample, they find that the stars of component 
C are a-enhanced, while the metal-rich stars of component A have 
[a/Fe] ~ 0. 1 , close to solar. They estimate isochrone ages for these 
stars, and find that the component C stars are almost all old, around 
10-12 Gyr, while the component A stars cover a wide age range 
from 2 to 14 Gyr. Their ages are consistent with the identifications 
of component A with the thin disk in a boxy/peanut bar structure, 
which has continued to undergo star formation, and with C being 
the old thick disk. Bensby et al, however, find only a small frac- 
tion of stars in the [Fe/H] range of what we believe to be the main 
bar/bulge (component B), around [Fe/H] = —0.25. 

We note that the iBensbv etaflfcoill) stars are mostly closer 
to the plane than our lowest fields at b = —5°. At b = —5°, at 
a radial distance of ± 0.5 kpc from the centre of the bulge, about 
13% of our stars have [Fe/H] between —0.2 and 0, but there is only 
one microlensed dwarf out of 26 in this abundance range. Along 
the minor axis at b = —5°, 40% of our stars have [Fe/H] > 0, 
somewhat lower than the 50% of microlensed stars in the same 
abundance interval which lie closer to the plane. This increase in 
the fraction of metal-rich stars closer to the plane is consistent with 
our finding that the fraction of the most metal-rich component A 
increases towards the plane. The number of stars in the metallicity 
range of component B (0 > [Fe/H] > —0.5) increases by a factor 
of only 1.2 from b — —10° to b = —5°, whereas the the number 
of stars in the metallicity range of component A ([Fe/H] > 0) in- 
creases by a factor of 5.6. Furthermore, if A is associated with the 
thin disk, then there may be an increasingly mixed-age population 
of metal-rich stars at low latitudes that is not seen at higher latitudes 
(b < -5°). 

The mi s sing p opulation B at the latitudes measured by 
IBensbv etailfcO 111) suggests that the component B stars close to 
the plane are greatly outnumbered by the metal-rich stars of com- 
ponent A. The region near the plane may be dominated by the 
flatter bar and host a younger, thin disk population which may be 
part of the bar. The low fraction of stars from component B in the 
IBensbv et al] J201lh data suggests that the stars of the early thin 
disk with [Fe/H] « —0.25 were mostly displaced to z— heights 
> 0.7 kpc during the boxy/peanut-bulge buckling process. As a 
check on these speculations, it would be interesting to know the 
abundance distribution of microlensed dwarfs at higher latitudes, 
closer to ours. The probability of finding microlensed stars de- 
creases rapid l y, how ever, with latitude. It is also possible that the 
IBensbv et af] ( 120 111) stars do not lie in the bulge itself: the geom- 
etry of microlensing would favor their preferential location on the 
far side of the peanut/boxy bulge. In this case, one would expect to 
find relatively few stars of component B (the dominant peanut/boxy 
bulge population) in the microlensed sample. To check this sugges- 
tion, we would need a realistic calculation of the dependence of 
microlensing optical depth on the position of stars in the inner re- 
gions of the Galaxy. 
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Component 


ARGOS 2012 




Bensby 


Babusiaux 


Hill et 








et al. 


et al. 


al. 2011 








2010 


2010 




b -5° 


-7.5° 


-10° 


< -5° 


-6° 


-4° 


A +0.12 


+0.11 


+0.08 


+0.32 


+0.14 


+0.32 


B -0.26 


-0.28 


-0.30 




-0.27 


-0.30 


C -0.66 


-0.68 


-0.70 


-0.60 






D -1.16 


-1.20 


-1.19 




-1.09 




E -1.73 


-1.67 


-1.68 









Table 6. The components of [Fe/H] identified in ARGOS survey compared 
to Bensby, Babusiaux and Hill surveys. The ARGOS survey values are 
taken for the integrated MDF's, across (=±15°, for Rq < 3.5 kpc. The 
metallicity of A is 0.15/0.13/0.09 at b = -5°, -7.5° and b = -10°, tak- 
ing only minor axis fields. 



9 CONCLUSIONS 

By measuring stellar parameters for the 25,500 stars in our survey, 
we are able to identify the nearby dwarfs and subgiants, determine 
distances to the giants and make a distance cut to separate out ap- 
proximately 14, 150 bulge candidates with \Rg\ ^ 3.5 kpc. The 
metallicity distribution function of these 14, 150 stars shows four 
to five distinct components (denoted A to E). We find that the rela- 
tive contributions of these components changes with position in the 
bulge. Our data extend into the disk beyond the bulge, which helps 
to identify the nature of the MDF components and interpret how 
they originated from the early disk of the Milky Way. 

The components show that the stars of the inner Galaxy are 
a composite population. The inner region of the Milky Way com- 
prises stars with [Fe/H] > —0.5, which appear to belong to the 
boxy/peanut bulge, and stars with [Fe/H] < —0.5 which lie in the 
inner regions but are not part of the boxy/peanut structure. There 
are two spatially and chemically separated bulge components with 
[Fe/H] > —0.5: (i) component A: a relatively thin boxy/peanut- 
bulge with mean [Fe/H] « 0.15 which we associate with the thin 
disk and (ii) component B: a thicker boxy/peanut-bulge population 
with mean [Fe/H] w -0.25. 

The duality of the bulge is clear but the mechanism by which 
this duality has arisen is not so clear. N-body models of disk evolu- 
tion show that the thick boxy/peanut component B is likely to have 
formed in a major instability event which occurred many Gyr ago. 
The thinner boxy/peanut bulge component A appears also to be as- 
sociated with such an instability event. It may have originated in 
the same event from the initial stratification of the thin disk, or it 
might have formed after the thicker structure. 

Chemically we associate the thin boxy/peanut-bulge (com- 
ponent A) with the younger thin disk, and this component can 
be followed out to large Galactic radii beyond the bulge at low 
| z\— heights from the plane where it merges with component B. 
It shows a radial metallicity gradient, becoming more metal-poor 
in the outer regions. Its gradient on the near side of the bulge from 
Rq = to 6 kpc is —0.02 ± 0.01 dex/kpc. Spatially it seems likely 
that component A comes from the redistribution of the younger 
thin disk in the inner regions into a boxy/peanut structure, via the 
bar-forming and bar-buckling instabilities of the disk. Our MDFs 
indicate that nearer the Galactic plane this population will be very 
dominant. It may also be mixed in age, due to ongoing star forma- 
tion, with younger stars located nearest to the plane. 

Component B is the primary bulge component that is more 
vertically extended than the metal rich component A. Outside the 
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bulge at high latitudes, where very few metal rich stars are present 
even inside the bulge, it represents the most metal rich fraction of 
stars observed in our sample at large |.z| -heights from the plane. It 
is associated with the early disk and from the instability scenario 
represents the primary population which formed the boxy/peanut 
bulge. 

The distinct mean [Fe/H] and [a/Fe] values for components 
A and B suggest that they formed in two distinct star forming 
episodes. Comparison of vertex deviations observed in the inner 
Galaxy and derived for the N-body models indicate that the thin 
boxy/peanut bar component A and the thick boxy/peanut bulge 
component B may lie at different angles to the Sun-center line: 
45° and 20° respectively. The two components A and B show the 
same vertical abundance gradient of: —0.08 ± 0.05 dex/kpc and 
—0.08 ± 0.04 dex/kpc, respectively. An even smaller vertical abun- 
dance gradient of —0.04 ± 0.02 dex/kpc is seen for the thick disk 
component C in the inner Galaxy. These small vertical gradients 
need not be interpreted as favouring a merger scenario for the for- 
mation of the bulge, or formation via dissipational collapse. They 
are likely to be residual gradients of the disk surviving from the in- 
stability event. The survival of vertical gradients in the initial disk 
needs to be confirmed by a more detailed study of the way in which 
the instability process in the N-body models maps the initial phase 
space into the post-instability phase space. We argue that the much 
larger overall vertical abundance gradient in the bulge region, of 
about —0.45 to —0.6 dex kpc -1 , comes primarily from changing 
contributions of components A-C with height above the plane (see 
Figure fTTTl. 

We find that the bulge is not a single population. Its metallic- 
ity components reflect the Galactic components which were present 
before the boxy/peanut bulge formation, and the properties of these 
components reflect the dynamical processes that altered the spatial 
distribution of stars primarily in the early disk. Stars with [Fe/H] 
> —0.5 are involved in the boxy/peanut bulge structure, and stars 
with [Fe/H] < —0.5 are present in the inner bulge regions and con- 
tribute to the total population of the inner Galaxy, particularly at 
higher latitudes. Our thin and thick boxy/peanut-bulge/bar struc- 
tures are both similar to the boxy/peanut bar/bulge seen in our in- 
stability driven N-body model, but with different scale heights and 
correspondingly different kinematics. Our more metal-poor com- 
ponents C, D and E, with [Fe/H] < —0.5 can be identified with the 
pre-instability thick disk, metal-weak thick disk and halo popula- 
tions. 

The presence of the two kinematically and chemically distinct 
components A and B in the boxy/peanut bulge is evidence that in- 
ternal evolution is the dominant process in the formation of the 
bulge. It seems unlikely that major mergers have been an impor- 
tant part of the formation of the bulge of the Milky Way. Such dis- 
ruption to the inner region would probably smooth the two bulge 
components into a continuous population. The survival of the two 
components of the bulge constrains the level of any major distur- 
bance that could have occurred in the inner Galaxy after the disk 
instability event. At this stage, we do not attempt to place limits on 
the contribution of any merger-generated "classical" bulge compo- 
nent. Further N-body modelling is needed to quantify this issue. 
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10 APPENDIX 

Regarding our red clump magnitude calibration: the Alves (2000) 
calibration is for stars in the range 0.0 > [Fe/H] > 0.5, which is the 
metallicity range of the majority of our stars in the bulge. Alves in- 
vestigated the metallicity dependence on the absolute magnitude of 
the clump and recommended n o [Fe/H] correction be made for the 
stars in this abundance range. lLanev. Joner. & PietrzvnskH d2012l) 
similarly reported no dependence on metallicity and an abso l ute re d 
clump magnitude of Mk = -1.61, in agreement wifh lAlvesI d2000l) . 

We have adopted his Mk over a broader metallicity range, 
because there is disagreement in the literature about the existence 
of a K-band magnitude-metallicity dependence. It is consistently 
reported to be weak at most, and there is disagreement about the 
factor in the equation for the metallicity dependence. We believe 
that it is appropriate to adopt a me an value of -1.61 for calcu lating 
clump distances for our stars (see IValentini & Munarill201fl for a 
review). Even if a metallicity dependence exists, its effect on our 
distances is likely to be small. For example, with Udalski's I-band 
calibration of a 0. 1 3 mag dex~ 1 dependence on abundance, the dis- 
tance of a star with [Fe/H] = -1.0 would change by about 6%. We 
note that 95% of our stars have [Fe/H] > - 1.0. 

Pietrzyriski, Gicrcn. & Udalski (2003) looked at the red clump 
magnitude dependence on age and [Fe/H] in J and K and found only 
a low or no dependence of K magnitude on [Fe/H] over the wide 
metallicity range of the LMC, SMC and the Carina and Fornax 
dwarf galaxies. They concluded that the mean K band magnitude 
of the red clump is an a ppropriate distance indicator for a wide 
range in [Fe/H] and age. Ivan Helshoecht & Groenewegenl d2007l) 
similarly found that the absolute magnitude is an appropriate dis- 
tance indicator in the range 0.4 > [Fe/H] > 0.5, although they 
adopted a slightly fainter value of Mk = —1.57 for the mean ab- 
solute magnitude of the clump. 

We note that theory predicts a metallicity dependence o f the 
absolute magnitude of the clump (i.e. lGirardi & S alaris 200 j|) and 
this has been supported by earlier observations (i.e. Saraiedini 
1 19991 : iTwarog. Anthonv-Twarog. & Brickej[l999l) . lUdalskil Jl99^) 
found a weak metallicity dependence in the I-band in the range 
0.6 < [Fe/H] < 0.2 using red clump stars in 15 Magellanic cloud 
clusters. However, he concluded that the absolute I-band magni- 
tude of the red clump is independent of age for ages 2-10 Gyrs. In 
summary, we will adopt Mk = —1.61 for the clump, independent 
of age and abundance. 
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